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ABSTRACT 



We present the results of a study of weak gravitational lensing by galaxies using imaging data that were obtained as 
part of the second Red Sequence Cluster Survey (RCS2). In order to compare to the baryonic properties of the lenses 
we focus here on the ~300 square degrees that overlap with the data release 7 (DR7) of the Sloan Digital Sky Survey 
(SDSS). The depth and image quality of the RCS2 enables us to significantly improve upon earlier work for luminous 
galaxies at z > 0.3. To model the lensing signal we employ a halo model which accounts for the clustering of the lenses 
and distinguishes between satellite and central galaxies. Comparison with dynamical masses from the SDSS shows 
a good correlation with the lensing mass for early-type galaxies. The correlation is less clear for late-type galaxies, 
possibly due to rotation. For low luminosity (stellar mass) early-type galaxies we find a satellite fraction of ~40% 
which rapidly decreases to < 10% with increasing luminosity (stellar mass). The satellite fraction of the late- types has 
a value in the range 0-15%, independent of luminosity or stellar mass. At high masses the satellite fraction is not well 
constrained, which we partly attribute to the modelling assumptions. To infer virial masses we apply simple models 
based on an independent satellite kinematics analysis to account for intrinsic scatter in the scaling relations. We find 
that early-types in the range 10^'' < Lr < 10^^'^Lq have virial masses that are about five times higher than those of 

2 3.+O.09 

late-type galaxies and that the mass scales as M200 oc L ' -oie, Por an early-type galaxy with a fiducal luminosity of 
lO^^I/r,0, we obtain a mass Af200 = (1.931q 14) x lO^^h^^ Mq. We also measure the virial mass-to-light ratio, and find 
for 1/200 < 10^^ Lq a value of M2(io / L200 = 42± 10 for early-types, which increases for higher luminosities to values that 
are consistent with those observed for groups and clusters of galaxies. For late- type galaxies we find a lower value of 
Af200 7^200 = 17 ± 9. Our measurements also show that early- and late-type galaxies have comparable halo masses for 
stellar masses Af* < 1Q" Mq, whereas the virial masses of early-type galaxies are higher for higher stellar masses. To 
compare the efficiency with which baryons have been converted into stars, we determine the total stellar mass within 
r-200. Our results for early- type galaxies suggest a variation in efficiency with a minimum of ^^10% for a stellar mass 



Af. 



,200 



— 10 Mq. The results for the late-type galaxies are not well constrained, but do suggest a larger value. 
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1. INTRODUCTION 

There is now overwhelming evidence that galaxies are 
surrounded by dark matter haloes. Studying the global 
properties of the haloes, such as their virial masses or den- 
sity profiles, however, has proven difficult due to a lack 
of reliable tracers of the gravitational potential at large 
distances. Improving observational constraints is impor- 
tant because the details of galaxy formation are not com- 
pletely clear, even though significant progress has been 
made in recent years (e.g. iBower et all 120101 : iKim et aLl 
120091) . The relation between the baryons and the dark mat- 
ter i n galaxies has been studied using numerical simulations 
'e.g.lWang et al.ll 2006: Cr oton et all2006l: ISomerville et al.l 



l2008t iMoster etal.. .2010:) and it is important to confront 
the predictions with observations. This requires reliable es- 
timates of both the dark matter and the baryonic content 



of galaxies. 

Several observables can be used to trace the baryons, 
such as the luminosity of a galaxy, which is readily avail- 
able. It is also possible to derive stellar masses by fitting 
stellar sy nthesis models to either the spectral feature s of 
a galaxy (iKauffmann et ahl 20031: iGallazzi et al.l 20051) or 
to its colours (|Ben fc de Jong|l200lHSalim et al.ll2007D The 
stellar mass estimates are tightly correlated to various other 
important global properties of galaxies (colour, mctallicity, 
luminosity, environment, see e.g. [Griitzbauch ct al. 20ll|, 
and references therein) and they are therefore considered a 
useful tracer of the baryonic content of a galaxy. 

Numerical simulations suggest that the dark matter 
haloes of massive gal axies extend out to hundreds of kilo- 
parsecs (e.g. Springel et al.l 2003). wh ich is supported by 
observations (e.g. iHoekstra et 311120041) . For nearby galax- 
ies it is possible to study the dark matter distribution using 
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the d ynamics of planetary nebulae (e.g. iNapolitano et al] 
l2009t) . In addition, studies of satellite galaxies around cen - 
tral galaxies (e.g. iMore et all 120111: IConroy et all |2007[ ) 
have provided constraints on the relation between baryons 
and dark matter. Unfortunately these studies require spec- 
troscopy of large numbers of objects, which makes them 
rather expensive. Furthermore, the observations are lim- 
ited to small scales due to the requirement of having optical 
tracers, which complicates the determination of the virial 
mass of the haloes galaxies reside in, unless one is willing 
to extrapolate the measurements. 

Fortunately it is possible to probe the matter distri- 
bution on large scales, thanks to an effect called weak 
gravitational lensing; we can measure the distortion of the 
shapes of faint background galaxies (sources) caused by the 
bending of light rays by intervening mass concentrations 
(lenses). The distortion is independent of the type of mat- 
ter in the lenses, and so the projected mass of the lens is 
measured without any assumption on the physical state of 
the matter at scales from a few kiloparsec to a few mega- 
parsec. 

The weak lensing signal around a single galaxy is too 
weak to detect since it is 10-100 times smaller than the in- 
trinsic ellipticities of galaxies. Therefore the galaxy-galaxy 
signal has to be averaged over many lenses to decrease 
the shape noise. Although individual galaxies cannot be 
studied in this way, their ave r age p r operties can be de- 
termined (e.g. iBrainerd et all 1 19961 iFischer et aD l2000t 
iHoekstra et al.ir2004 ). Only more recently has it become 
possible to study lenses as a function of properties such 
as type, luminosity, stellar mass, etc., because early studies 
lacked the ancillary data needed to subdivide th e lenses into 
subsamples. For instance lHoekstra et al.l ()2005l ) used nearly 
34 squ are degrees of the Re d Sequence Cluster Survey 
(RCS) (iGladders fc Y cc' 2005? ) for which photometric red- 
shifts were available (jHsich et al.l 120051 ). to study the re- 
lation between the virial mass and baryonic contents of 
isolated galaxies in the redshift range 0.2 < z < 0.4, 
and derived star formation efficiencies for early- and late- 
type galaxies. Thanks to the wealt h of ancillary data , the 
Sloan Digital Sky Survey (SDSS; lYork et al.l l2000h has 
had a major irnpact on galaxy-galaxy lensing studies (e.g. 
iGuzik fc Selj al? 2002; 'Mandclbaun ]_et al.l[2006h . This is ev- 
idenced bv lMandclbaum ct al. ( 20061) who used nearly 5000 
squar e degrees of the SDSS DR4 (lAdelman-McCarthv et al.l 
l2006t ) to study galaxies in the redshift range 0.02 < z < 0.35 
as a function of galaxy type and environment, and con- 
strained the stellar mass to virial mass relation, the lumi- 
nosity to virial mass relation and the satellite fractions of 
the lens samples. 

Currently no survey can surpass the precision that can 
be achieved by the SDSS at low redshift {z < 0.3) because 
of the large survey area and the availability of spectroscopic 
data. We note, however, that complementing the SDSS data 
with deeper imaging by the Panoramic Surv ey Telescope 
fc Ra pid Response Systen:(3 (Pan-STARRS; iKaiser et aD 
|2002|) will provide a major improvement, as is demonstrated 
by the results we present here. For lenses with z > 0.3 it 
is possible to achieve a significant improvement over the 
SDSS results by surveying a smaller area with deeper data 
and good image quality; it allows us to use sources at higher 
rcdshifts. This is important because the amplitude of the 



lensing signal scales proportionally to the ratio of the an- 
gular diameter distance between the lens and the source 
and the distance between the observer and the source. The 
signal decreases rapidly when the lens redshift approaches 
the peak of the source redshift distribution, which occurs 
around z - 0.35 for the SDSS. 

In this paper we use data from the second generation 
Red Sequence Cluster Survey (RCS2) to measure the weak 
lensing signal around galaxies that are observed in the 
SDSS. The RCS2 is a nearly 900 square degree imaging 
survey carried out by the Canada-Francc-Hawaii- Telescope 
(CFHT), and is ~2 magnitudes deeper than the SDSS in 
r' . The increase in depth combined with a median seeing 
of 0.7", which is a factor of two smaller than the seeing in 
the SDSS, results in a source galaxy number density that is 
about five times higher, and a source redshift distribution 
that peaks at Z'-^^O.?. 

We use the overlapping area between the two surveys, 
which amounts to approximately 300 square degrees, in or- 
der to assign the spectroscopic redshifts, luminosities, stel- 
lar masses and dynamical masses from the SDSS to the 
lenses. The lensing analysis itself is performed on the RCS2 
data. Even though the overlap between the surveys is mod- 
est, the loss in survey area is outweighed by the gain in the 
number density of source galaxies and the improvement of 
the lensing efficiency. This enables us to improve the mea- 
surements of the lensing signal around the most massive 
galaxies, which mostly reside at redshifts where the SDSS 
is not very sensitive. 

In this paper we describe the lenses in Section [51 The 
weak lensing analysis is discussed in Section [31 The halo 
model that we have implemented is introduced in Section 
[31 In Section[Slwe compare the weak lensing mass to the dy- 
namical mass. We describe the luminosity results in Section 
[6l and the stellar mass results in Section [71 Wc summa- 
rize our conclusions in Section [SI Throughout the paper 
we assume a WMAP5 cosmology (jKomatsu et al.l l2009t) 
with erg = 0.8, flA = 0.73, nM = 0.27, f^b = 0.045 and 
the dimensionless Hubble parameter h = 0.7. All distances 
quoted are in physical (rather than comoving) units unless 
explicitly stated otherwise. 



2. LENS SAMPLE 

The SDSS has imaged roughly a quarter of the en- 
tire sky, and has observed the spect r a for about one mi l- 
lion galaxies (jEisenstein et al.l I2OOII : IStrauss et al.ll2002| ). 
The combination of spectroscopic coverage and photome- 
try in five optical bands {u,g,r,i, z) in the SDSS provides 
a wealth of galaxy information that is not available for the 
RCS2. To use this information, but also benefit from the im- 
proved lensing quality of the RCS2, we use the 300 square 
degrees overlap between the surveys for our analysis. We 
matc h the RCS2 catalogues to the DR7 (jAbazaiian et alJ 
I20090 spectroscopic catalogue, to the MPA-JHU DR713 
stellar mass catalogue and to the NYU Value Added 
Galaxy Catalogue (NYU-VAGCf[ (iBlanton et al. 



2005 



lAdelman-McCarthv et~alll2008l: iPadmanabhan et al.ll200 
which yields the spectroscopic redshift, luminosity, stellar 
mass, and the dynamical mass of 1.7 x 10^ galaxies. These 
form the lens sample of this work; we study the distortion 



http://pan-starrs.ifa.hawaii.edu/public/ 



^ http://www.nipa-garching.mpg.de/SDSS/DR7/ 
'' http://sdss.physics.nyu.edu/vagc/ 
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these galaxies imprint as a function of their baryonic con- 
tent on the shapes of the background galaxies. 

As the relation between dark matter and baryons de- 
pends on galaxy type, we split the lens sample into early- 
and late-type galaxies using the frac-deV parameter in- 
cluded in the SDSS photometric catalogues. This parame- 
ter is determined by simultaneously fitting frac-deV times 
the best-fitting De Vaucouleur profile plus {l-frac-deV) 
times the best-fitting exponential profile to an object's 
brightness profile. This has been done in the g, r and i 
band, and we use the average value. We classify galax- 
ies with frac-deV > 0.5 as early-types, and galaxies with 
frac-deV < 0.5 as late-types. The classification of early- 
types is at least 96% complete and 76% reliable (96% of 
all early- type galaxies are in the early- type sample, while 
76% of all the galaxies in the early-type sample are actu- 
ally early-types), and the classificati on of late-types is at 
least 55% complete and 90% reliable (jStrateva et al.ll200lt 
iMandelbaum et al.|[2006t) . 

We visually inspect the brightest and most massive 
early- and late-type galaxies of our lens sample using our 
RCS2 imaging data. We find that about 30 of the 100 
most massive late-types (with a stellar mass in the range 
XO^i-^— IO^^-^Mq) actually consists of multiple objects with 
small separations. These galaxies reside at a redshift of 
~ 0.4, and are not well resolved in the SDSS. They are 
not removed from the analysis as that may introduce a se- 
lection bias. More importantly, including them facilitates 
a comparison to the literature. As a test, we excluded 
these lenses, and found that the results did not significantly 
change (note, however, that due to the low number of mas- 
sive late-type lenses, the errors are large). 

2.1. Luminosities & Stellar Masses 

The MPA-JHU stellar mass catalogue contains about 
7 X 10^ unique galaxies, and provides the r-band absolute 
magnitudes and the stellar mass estimates of our lenses. 
The absolute magnitudes that are used to compute the lu- 
minosities and stellar masses are based on the Petrosian 
apparent magnitudes from the SDSS. The Petrosian appar- 
ent magnitude measures the fiux within a circular aperture 
whose radius depends on the azimuthally averaged bright- 
ness profile in the r-band. It does not include the flux at 
very large radii from a galaxy, and therefore underestimates 
the total flu x by typically a few tenths of a magnitude 
(JBlanton et a l. 2001). Although we do not correct for the 
missing flux as it would complicate a comparison with pre- 
vious observational work, this should be kept in mind when 
comparing our results to predictions from numerical simu- 
lations. 

The absolute magnitudes have been corrected for ex - 
tinction using the dust maps from ISchleeel et al.l (jl998t) , 
the k-corrections have been calc ulated to z = 0.0 
using the KCDRRECT v4_2 code (jBlanton et al.l l2003t 
iBlanton &: Roweisll2007() . and the distance modulus is de- 
termined with h = 0.7. We convert the absolute magnitudes 
into solar luminosities using the absolute AB magnitude in 
the SDDS r-band of Msoiar = 4.65 for z = 0.0. We ac- 
count for passive evolution by dividing the luminosities of 
the early-type galaxies by {1 + z). The luminosity evolution 
of late-type galaxies can in principle be computed if the 
star formation histories (SFHs) are accurately known. The 
SFHs are generally uncertain, however, since they depend 
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Fig. 1. Stellar mass versus luminosity of the lens sample. 
The colour coding represents the redshift of the galaxies 
as denoted in the lower right-hand corner. The histograms 
for all galaxies (black line), the early-types (red line) and 
the late- types (blue line), as a function of stellar mass and 
luminosity are also shown, and are drawn slightly offset 
for clarity. The dashed diagonal lines indicate the addi- 
tional mass-to-light ratio cuts we have applied (objects with 
Mtf/Lr between 0.2 and 10 have been selected) to remove 
outliers that may contaminate the lensing signal. 



on many parameters such as the stellar mass, environment, 
assembly history, and AGN activity of a galaxy. Hence the 
luminosity evolution is difficult to determine and the cor- 
rection highly uncertain. We therefore do not correct the 
luminosites of late-type galaxies for evolution. 

The stellar masses have been estimated by fitting a li- 
brary of iBruzual fc CharlotI (l2003l ) stellar population mod- 
els to the u,g,r,i,z photometry of the galaxies in the 
S DSS. The initi al mass function (IMF) was taken to be 
a iKrouDa' (2001) IMF and the modelling methodology fol- 
lows Sahm et al. (2007[ ) . 

Nearly all galaxies with a spectroscopic redshift from 
DR7 are present in the stellar mass catalogue. Figure [1] 
shows the stellar mass versus luminosity for the matched 
galaxies. The different colours represent galaxies at differ- 
ent redshifts. The most massive and luminous galaxies in 
our sample reside in the highest redshift range, and are 
almost exclusively early-type galaxies. Also shown are the 
histograms of the stellar masses and of the luminosities on 
respectively the horizontal and vertical axis. The dashed 
lines indicate the additional 0.2 < M^/L^ < 10 cut we ap- 
ply to minimize the outlier contamination of the lensing 
bins. 
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2.2. Dynamical Masses 

The motions of stars in a galaxy provide an alternative 
way to estimate the mass of a galaxy at small radii, and 
constrain the scaling relations between baryons and dark 
matter. Spectroscopic observations are required to measure 
the velocity dispersion, which is converted into a dynamical 
mass estimate via the scalar virial theorem, taking into ac- 
count projection effects and assumptions on the structure 
of the stellar orbits: 



GAfdyn = Kvin)af^^Re 



(1) 



with (Tios the line-of-sight velocity dispersion of the galaxy. 
Re the effective radius (containing 50% of the light of the 
best fit Sersic model), and Kv{n) a term that includes the 
effects of struct ure on stellar dyn amics, which can be ap- 
proximated by (jBertin et al.ll2002n : 



Kvin) 



73.32 



10.465+ (71-0.94)2 



0.954, 



(2) 



with n the Sersic index ()Sersiclll968[ ). 

Using the dynamical mass as a tracer for the total 
mass of a galaxy has various complications. Firstly, it is 
implicitly assumed that the velocity dispersion in Equation 
[T]is only generated by the radial motions of the stars, and 
the Kvin) term is derived under the assumption that the 
mass distribution is spherical, dynamically isotropic, and 
non-rotating. In reality, however, the rotation of a galaxy 
contributes to the measured velocity dispersion as well, and 
this effect is particularly important in late-type galaxies. 
The majority of the early-type galaxies in our study are 
mass ive and luminous. Th ey are expected to rotate slowly 
(e.g. lEmsellem et al.ll2007t ). so their dynamical mass esti- 
mates are less affected. The dynamical masses of late-type 
galaxies, however, are potentially overestimated. A second 
complication arises from the fact that the spectroscopic 
fibre within which the velocity dispersion is measured has 
a fixed size. Therefore, the physical region over which the 
velocity dispersion is averaged depends on the redshift of a 
galaxy, and hence it probes different regions for galaxies at 
different rcdshifts. If the velocity dispersion changes with 
radius, we would effectively assign different dynamical 
masses to the same galaxy depending on its redshift. 
Thirdly, the dynamical mass is measured within the 
effective radius. The effective radius is a rather arbitrary 
point, as it depends on parameters such as the shape, the 
brightness profile and the orientation of a galaxy, and the 
distribution of dust within the galaxy. Even if a galaxy is 
spherical and isotropic, it is not clear whether the effective 
radius marks a special point in relation to the total mass 
content of a galaxy, given that the dark matter does not 
follow the distribution of stars. This is most obvious in 
the outer regions of a galaxy, where most of the matter is 
dark. 

To calculate the dynamical mass of our lenses, we re- 
trieve the velocity dispersions from the SDSS spectroscopic 
catalogue. As it is complex to estimate the velocity disper- 
sion of galaxies whose spectra are dominated by multiple 
components, e.g. galaxies with different stellar popula- 
tions or different kinematic components, the SDSS only 
provides estimates for spheroidal systems whose spectra 
are dominated by red stars. At low redshift, the selection 
also includes the bulges of late-type galaxies because their 



spectra are similar to the spectra of early-type galaxies. 
The Sersic index and the effective radius are obtained from 
the NYU-VAGC. The sizes and fluxes are underestimated 
10% and 15% respectively for large galaxies an d galaxies 
with high Sersic indices (iBlanton et al.l l2005t ). whereas 
the Sersic index itself is underestimated by ^ 0.5 to 
~ 1.3 for g alaxie s with high Sersic indices. It is shown in 
iGuo et ahl (J2009t) that these biases arise from background 
overestimation and subtraction. As a result, the dynamical 
mass estimates for these galaxies may be slightly biased, 
but we do not account for it since we do not know the 
correction for each galaxy. To ensure that the dynamical 
mass is computed in approximately the rest-frame r-band, 
we split the sample according to redshift. For galaxies at 
z < 0.2 we use the Sersic index and effective radius in 
the r-band, for galaxies between 0.2 < z < 0.4 we use the 
values in the i-band, and for galaxies at z > 0.4 we average 
the values of the i- and z-band. 



3. LENSING ANALYSIS 

3.1. The RCS2 

The lensing signal can be detected with high signifi- 
cance at low redshifts (z < 0.3) using SDSS data only. At 
higher redshifts, the significance decreases rapidly, because 
of the limited imaging depth and image quality of the SDSS. 
To improve the lensing signal-to-noise ratio at z > 0.3, 
we use the deep imagin g data from t he Red Sequence 
Cluster Survey 2 (RCS2) (icilbank et al.l[20Tl[) instead. The 
RCS2 is a nearly 900 square degree imaging survey in three 
bands (g', r' and z') carried out with the Canada- France- 
Hawaii Telescope (CFHT) using the 1 square degree cam- 
era MegaCam. The primary survey area is divided into 13 
well-separated patches on the sky (including the uncom- 
pleted patch 1303), each with an area ranging from 20 to 
100 square degrees Q Since the RCS2 consists of single ex- 
posures only, it is difficult to identify cosmic rays, especially 
those that hit stars and galaxies. However, only a small 
fraction of objects is hit by a cosmic ray, and the affected 
objects do not bias the measurements, bu t act as a negligi- 
ble source of noise (jHoekstra et al.ir2004|) . We perform the 
weak lensing analysis in the SDSS and RCS2 overlap using 
the 8 minute exposures of the r-band {r'^^ ~24.8), which 
is best suited for lensing as it has a median seeing of 0.7". 



3.2. Image processing 

We retrieve the ElixiiQ processed images from the 
Canadian Astronomy Da ta Centre (CADC) archiv Q We 
use the THELI pipeline (JErben et al.ll2005ll2009t) to sub- 
tract the image backgrounds, to create weight maps that 
we use in the object detection phase, and to identify satel- 
lite and ast eroid trails. T o obtain accurate astrometry, we 
run SCAMP (jBertinl [20061 ) on the images, which enables us 
to match our catalogues to the SDSS. The polynomial co- 
efficients from SCAMP describing the mapping from image 



* The CFHT Legacy Survey Wide, comprising of 171 square 
degrees of imaging data in m*, g', r', i' and z', is also included 
in the RCS2, but is not used in this study. 
^ http://www.cfht.hawaii.edu/Instruments/Ehxir/ 
" http://wwwl.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/cadc/ 
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Flux radius [pix] 

Fig. 2. The size- magnitude diagram of one of the chips 
in a randomly picked exposure. The black dots are the 
SExtractor detections, the green dots are the selected 
stars, the red dots are the 3-sigma outliers, and the blue 
dots are the selected galaxies. The dashed lines indicate 
the location of the stellar branch. Thanks to the good im- 
age quality the stars are easily separated from the galaxies. 



to sky coordinates are used to calculate the camera dis- 
tortion. We use the automated masking routines from the 
THELI pipeline to generate image masks and to combine 
them with the RCS2 masks in order to omit image regions 
that contaminate the lensing signal (e.g. saturated stars, 
satellite trails). All masks are inspected by eye, and manu- 
ally improved where nec essary. 

Wc use SExtractor (JBertin fc Arnoutsl[T996[) to detect 
the objects in the images. To select the stars for modelling 
the PSF variation across the images, we first identify the 
locus of the stellar branch in a size-magnitude diagram. We 
select the non-saturated objects close to the stellar branch 
with a signal to noise ratio larger than 30 and with no 
SExtractor flags raised. To remove small galaxies that have 
been misidcntified as stars, and stars that have been af- 
fected by cosmic rays, we fit a second-order polynomial to 
both the size and the ellipticity of these star-candidates, 
and discard all 3-sigma outliers. We clean the stellar se- 
lection even further in the shape measurement pipeline by 
removing shape parameter outliers. All objects larger than 
1.2 times the local size of the PSF are classified as galaxies. 

In Figure [2] we illustrate the star-galaxy separation. It 
has been fully automated, but as a precaution we inspect 
all size-magnitude diagrams by eye. The separation fails for 
a few chips that have either very few stars or a PSF with 
a large FWHM, and we manually adjust those. As neigh- 
bouring patches overlap by ~ 1 arcminute, we remove all 
galaxies within 35 arcseconds from the image edges in order 
to avoid duplicating the lenses and sources in our analysis. 
Elixir provides approximate zeropoints for each point- 
ing, which we use to measure the r'-band magnitudes of 




100 1000 

Projected separation [kpc] 

Fig. 3. The source galaxy overdcnsity as a function of dis- 
tance from the lenses for the different stellar mass bins. 
The overdensity increases with stellar mass. Massive galax- 
ies reside on average at higher redshifts and live in denser 
environments with more satellite galaxies. 



the objects in the images. We correct the magnitudes for 
galact ic extinction using the dust maps from lSchlegel et all 
(|1998[ ). These magnitudes are not a s accurately calibrated 
as those from iGilbank et al.l (J201l|) , and differ in the r'- 
band on average by —0.01 ± 0.32. Our calibration is, how- 
ever, sufficiently accurate to select the source galaxy sam- 
ple. For the calculation of the luminosity overdensity, which 
is discussed in Sect ion 16. 1[ we use the catalogues from 
IGilbank et al] ()20lH) instead. 



3.3. Contamination correction 

A fraction of the galaxies in the source catalogue is 
physically associated with the lenses. Since we lack red- 
shifts for the sources, we are unable to remove them. These 
objects are not lensed, and therefore dilute the lensing sig- 
nal. To estimate this contamination we measure fcg (r) , the 
excess source number density around the lenses. We show 
the overdensity around the lenses which have been divided 
into seven stellar mass bins (defined in Table [3]) as a func- 
tion of lens-source separation in Figure [31 The error bars 
are computed assuming that the number of source galaxies 
in each radial bin follows a Poisson distribution. The con- 
tamination increases with stellar mass, as massive galax- 
ies reside in denser environments and therefore have more 
satellite galaxies. Although the overdensity is shown inde- 
pendently of the lens galaxy type in Figure |31 we measure 
it for the early- and late-types separately in the science 
analysis presented in Section [SJ [5] and [T] Assuming that 
the satellite galaxies have random orientations, we correct 
for the contamination by boosting the lensing signal with 
a factor 1 + /cg(r). Note, however, that the contamination 
correction may be too small if satellite galaxies are pref- 
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erentially radially aligned in the direction of the lens. This 
type of intrinsic alignment has been studied with seemingly 
different results; some authors (e.g. lAgustsson fc BrainerdI 
l2006t iFaltenbacher et al]|2007t) who determined the galaxy 
orientation using the isophotal position ang les, have ob- 
served a stronger alignment tha n others (e.g. iHirata et al.l 
2004 Mandelbauni et ani2005a ) who used galaxy moments. 



Siverd et al.l 



2009[) and lHao et all (|201lh attribute the dis- 
crepancy to the different definitions of the position angle 
of a galaxy. As we measure the shapes of source galaxies 
using galaxy moments, we expect that intrinsic alignment 
only has a minor impact on the correction factor and hence 
can be safely ignored. 

Gravitational lenses do not only shear the images of the 
source galaxies, but also magnify the background sky. As a 
result, the flux of the sources is magnified, and the source 
galaxy number density is diluted. These combined effects 
are known as magnification bias, and it changes the source 
density around the lenses. The effect is negligible for the 
lensing study presented here. 

3.4. Shape measurement 

The measurement of the shapes of galaxies is central to 
any weak lensing analysis. The accuracy that is required 
depends on the science goal. For example, in cosmic shear 
studies aimed at constraining cosmological parameters, 
it is necessary to accurately correct the measured galaxy 
shapes for the anisotropic smearing of the PSF since the 
signal is small and very sensitive to any PSF residual 
systematic. In contrast, in the case of galaxy-galaxy lensing 
the signal is averaged over many lens-source pairs with 
random orientations, which removes most of the PSF 
systematics on small scales. 

For our lensing analysis we m easure the s hapes 
of galaxies with the KSB method ([Kaiser et all Il995t 
iLuppino fc Kaisen 119971 : iHoekstra et al. 11998), us i ng th e 



implementation described bv iHoekstra et aTl ( 



19981 1200011 . 



The measured galaxy shapes are corrected for smearing 
by the PSF under the assumption that the brightness 
distribution of stars can be described by an isotropic profile 
convolved with a small anisotropic kernel. Generally, the 
PSF is more complicated which may lead to biases. The 
version of KSB we use has been tested on simulated 
images as part of the Shear Testing Prograni r ne (ST EP) 
1 and 2 (the 'HH' method in iHevmans etall ()2006D and 
iMassev et al.l (|2007t) respectively). These tests have shown 
that the correction scheme works well for a variety of 
PSFs; in STEP2, the HH method underestimates the shear 
on average by 1-2% only. 

The mapping between the sky coordinates and the 
CCD pixels is slightly non-linear due to the camera optics, 
which causes an additional shear that needs to be cor- 
rected. We calculate the shear induced by this distortion 
using the polynomial coefficients from SCAMP describing the 
mapping from image to sky coordinates. The camera shear 
of MegaCam is shown in Figure H) The images of both the 
stars and the galaxies are sheared, with a value reaching 
1.5% at the corners of the images. At large lens-source 
separations, where the gravitational lensing signal is small, 
the camera shear dominat es the observed lensing signal. 
IHoekstra et al.l (|1998L l2000[ ) demonstrate that the observed 
shear is the sum of the gravitational shear and the camera 
shear. We therefore simply subtract the camera shear from 
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Fig. 4. Shear induced by camera distortion in the 
MegaCam imager. The camera shear is largest in the cor- 
ners of the mosaic, with values up to 1.5%. As the observed 
shear is the sum of the gravitational shear and the cam- 
era shear, we simply subtract the camera shear from the 
observed galaxy ellipticities to correct for it. 



the observed ellipticities of the galaxies to correct for it. 

To demonstrate the excellence of the RCS2 as a lensing 
survey, we measure the galaxy-mass cross-correlation func- 
tion in the exposures that significantly overlap with the 
SDSS (defined as having more than 30 matching objects). 
301 exposures of the total overlapping 350 meet this re- 
quirement, which after masking and exclusion of the image 
boundaries leads to an effective area of approximately 260 
square degrees. The galaxy-mass cross-correlation function 
measures the correlation between the galaxies and the 
surrounding distribution of (predominantly dark) matter. 
We compute it by measuring the azimuthally averaged 
tangential shear as a function of radial distance from the 
lens: 

(7*)W = ^, (3) 

where AI](r) = fi(< r) — S(r) is the difference between the 
mean projected surface density enclosed by r and the mean 
projected surface density in an annulus at r, and Egrit is 
the critical surface density 



^crit — 



D, 



47rG DiDu 



(4) 



with _D;, Ds and Dis the angular diameter distance to the 
lens, the source, and between the lens and the source re- 
spectively. 

Since we do not have redshifts for all galaxies we 
separate the lenses from the sou rces using magnitude 
cuts (see e.g. IHoekstra et al.ll2004[ ). Objects with 19.5 < 
nij.! < 21.5 are defined as lenses, and objects with 22.0 < 
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Fig. 5. The galaxy-mass cross-correlation function around 
7.3x 10^ apparent magnitude selected lenses measured with 
5.9x10^ sources. The black symbols are the tangential 
shear, the red symbols are the cross shear. The top axis 
shows the projected separation in physical units for the 
median lens redshift Zmed=0.34. The inset shows the sig- 
nal on a linear scale for small separations. The signal has 
been corrected for contributions from systematic shear, and 
boosted to account for source galaxy contamination. The 
dashed (dotted) line shows the best fit SIS (NFW), fitted 
to the shear on scales between 0.2 and 0.6 arcminutes. The 
clustering of galaxies causes excess shear at scales >1 ar- 
cminutes. 



rrir' < 24.0 arc sources. We discard objects with cUiptic- 
ities larger than 1, and objects that have a SExtractor 
flag raised. Using these selection criteria we find 7.3x10^ 
lenses and 5.9x10^ sources. The corresponding effective 
source number density is 6.3 arcmin"^, which is five times 
higher than the so urce density of 1.2 arcmin "^ used in 
the SDSS analysis (jMandelbaum et al.|[2"005al ). To obtain 
the approximate redshift distribution of the lenses and 
sources, we apply identical magnitude cuts to the photo- 
metric redshift catalogues of the Canada-France-Hawaii- 
Telescope Legacy Survey (CFHTLS) "Deep Survey" fields 
(jllbert et al.ll2006l ). We stack the signals of all the lenses in 
the RCS2, and azimuthally average them in radial bins. To 
remove the contributions of systematic shear (from, e.g., 
the image masks), we subtract the signal computed around 
random lenses from the signal around the real lenses. We 
measure the source galaxy overdensity as a function of lens- 
source separation, and boost the signal to correct for the 
contamination as outlined in Section 13.31 Figure [S] shows 
the tangential shear, and the inset shows the signal at small 
scales using a linear vertical scale. 

We also measure the cross shear around the lenses 
by rotating the background galaxies 45° and repeating the 
measurement. Gravitational lensing does not produce cross 
shear, and a non-zero signal indicates the presence of resid- 



ual systematics in the catalogues. We indicate the cross 
shear with the red symbols in the inset in Figure [5l and 
note that it is consistent with zero on all scales. 

For reference, we fit a singular isotherm al sphere (SIS) 
and a Navarro-Frcnk- White (NFW) profile (jNavarro et al.l 
1996J) to the tangential shear on scales between 0.2 and 
0.6 arcminutes (^60-180 kpc at the median lens redshift 
Zmed = 0.34). The SIS signal is given by 



lt,sis(r) = — = - 



1 



2r 



D, 2r' 



(5) 



where te is the Einstein radius and a the velocity disper- 
sion. We indicate the best fit SIS model with the dashed 
line in Figure \5\ The NFW density profile is given by 



p{r) 



ScPc 



(r/r,)(l + r/r,)2 



(6) 



with Sc the characteristic overdensity of the halo, pc the 
critical density for closure of the universe, and r^ = 
''200/cNFW the scale radius, with cnfw the concentration 
parameter. The NFW profile is specified by two free pa- 
rameters: the mass and the concentration parameter. Since 
numerical simulations have shown that the concentration 
depends on the mass and redshift of the halo, we can reduce 
the number of free parameters in the fit by adopting a mass- 
concentra tion relation . We u se the mass-concentration rela- 
tion from lDuffv et al.l ()2008[ ). which is based on numerical 
simulations using the best fit parameters of the WMAP5 
cosmology. It is given by 



Cnfw — 5.71 



M: 



200 



2 x 1012/i-iA'fp 



(1 + ^)" 



(7) 



with M200 the mass in units of h~^MQ. M200 is defined as 
the mass inside a sphere with radius r2oo, the radius where 
the density is 200 times the critical density pc- We use the 
median lens redshift Zmed = 0.34 for the stacked lenses 
in the NFW fit, and calculate the tangential shear profile 
using the analytical expression s prov ided by iBartelmannI 
(|1996D and I Wright fc BrainerdI (|2000D . The best fit NFW 
profile is indicated by the dotted line in Figure [5] 

It is clear that both the SIS and NFW profiles underes- 
timate the signal at scales larger than ^1 arcminute, which 
corresponds to ~300 kpc at the median lens redshift. The 
majority of galaxies live in clustered environments, and 
with gravitational lensing we measure the shear induced 
by neighbouring galaxy haloes as well. This excess lensing 
signal complicates a straightforward analysis of the data. 
The problem could be avoided by studying the lensing 
signal o n small scales around isolated galaxies (following 
Hockstra et al.l (|2005l) ). but this requires the availability 
of redshifts for all galaxies, which we do not have in the 
RCS2. Alternatively, the lensing signal can be modelled 
taking the clustering of the lenses into account, which 
enables the simultaneous study of the mass and of the 
clustering prop erties of the galaxies. This is inhe rent in 
the hafo model (lSeliakll2000l: iCoorav fc Shetbl[200l . which 
wc will use here. 

The lenses in a bin generally have a range of masses. 
The correct interpretation of the signal therefore requires 
knowledge of the distribution of the masses of the lens 
galaxies, an issue we return to at the end of Section 2] 
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4. HALO MODEL 

Galaxies form in the gravitational potential of dark 
matter haloes and therefore trace the large scale dis- 
tribution of matter in the universe. The quantity that 
describes the relation between galaxies and dark mat- 
ter is referred to as galaxy biasing. The description of 
galaxy biasing is non-trivial as the physics governing 
galaxy formation is complex, and the bias may depend on 
the dark matter halo mass, environment, scale and red- 
shift (e.g. ICresswell fc Percivall [20091: ICoupon eraI1l201lt 
iKovac et al.ll2011[ ). To gain insight into the relation between 
galaxies and dark matter the weak lensing signal around 
galaxies can be used, as it measures the correlation be- 
tween the galaxies and the surrounding dark matter distri- 
bution. These lensing measurements provide constraints for 
models of the large scale distribution of matter, which are 
commonly described with the power spectrum of the den- 
sity fl uctuations (e.g. iPeacock fc DoddslllQQo : ISmith et al.l 
120031) . For a given po wer spectrum, the len sing signal can 
be computed directly (jGuzik fc Seliakll200H) : 



X f dkkPik,x,O)J2{kr{x)0), 



(8) 



Da 



with X the radial distance (in a fiat universe, x = 
with a the scale factor and Da the angular diameter dis- 
tance), W^i(x) the normalized radial distribution of the 
lenses, /(x) = J^ dx'gix,x')W2ix'), with W2{x') the ra- 
dial distribution of the sources, and 



5(X,X') = 



D,D, 



l^ls 



Dsa{zL)' 



(9) 



P{k) is the power spectrum under consideration, and J2 
is the second Bessel function of the first kind. Instead of 
using a single power spectrum to describe the distribution 
of matter in the universe, it is beneficial to consider the 
various components that contribute, as is done in the halo 
model. This allows a simultaneous study of the halo masses 
of galaxies and of their clustering properties. 

In the halo model the mass distribution in the universe 
is described as a distinct number of dark matter haloes 
that are clustered. As the large scale spatial distribution 
of haloes is unlikely to affect the physics inside individual 
haloes, and vice versa, the description of the model can be 
separated into two steps: the halo mass function and the 
bias at large scales, and the halo occupation distribution at 
small scales. 

The large scale distribution of haloes can be described 
by the halo number density. In the Press-Schechter ap- 
proach (jPress fc Schechterlll974[ ) the dark matter haloes 
are assumed to form by spherical collapse. This, however, 
leads to a halo number density that overestimates the abun- 
dance of galaxies below the non-linear mass scale. Better 
agreement with numerical simulations of hierarchical struc- 
ture formation comes from t he assumption of ellipsoidal 
rather that spherical collapse (jSheth et al.ll2001l ) . The num- 
ber density of bound objects is generally written as 



nh{M,z)dM = ^^f{y)dv, 



(10) 



where nh{M, z) is the halo mass function which depends on 
the halo mass M and redshift z, and p is the mean matter 
density of the universe at redshift z. Unless explicitly stated 
otherwise we use M — A/200 ■ The peak height v is given by 



Sscjz) 

cr{M,z) 



(11) 



with 6sc{z) the critical overdensity required for spherical 
collapse at redshift z, and cr(M, z) the rms of the density 
fluctuation field on the scale R = (3M/47rp)^/'^, extrap- 
olated to z using linear t heory. In the case of ellipsoidal 
collapse, f{v) is given by (jSheth et al.ll200l[ ) 



fi,y)^Ail + iai^)-P)iy 



-l/2^-au/2 



(12) 



with a = 0.707, p = 0.3, and A ~ 0.13683 a constant 
that is determined by requiring J f{i')dh' — 1 (i.e. mass 
conservation) . 

How the haloes trace the mass is given by the halo- 
to-mass bias, which is defined as the ratio of the power 
spectrum of the halo distribution to the power spectrum of 
the matter distributi on. We use an ana lytical formula for 
the bias as given by ISheth et al.l (12001 ), but in corporate 
the adjustments described in [Tinker et al.l (|2005[ ): 



Kiy) = 1 + 



Vo-Ssc 



{avY 



(«!/)= + 6(1 -c)(l-c/2) J' 
(13) 



with a = 0.707, b = 0.35 and c == 0.80. The scale depen- 
dence of the bias is given by 



b'^{v^r) = h'^{v) 



1.49 



[1-H.17^„,(r) 
[l+0.69e™(r)]2-09' 



(14) 



where S,m{'r) is the matter correlation function, which 
in turn is the Fourier transform of t he no n-linear power 
spectrum Pnl(^) from ISmith et all (|2003f ). and r is the 
distance to the centre of the halo. 

To describe how the galaxies and dark matter are dis- 
tributed wi thin the haloes, we clo sely follow the approach 
outline d in iGuzik fc SeliakI (J2002t) and iMandelbaum et all 
(l2005b[ ). Galaxies living inside dark matter haloes are 
divided into two classes; they are either a central galaxy 
located in the central halo, or a satellite galaxy located in 
a subhalo inside the central halo. The fraction of satellites 
in a certain sample of galaxies is denoted by a. The 
number of satellites in a central halo is described by the 
halo occupation distribution (HOD). Galaxy formation 
simulations (e.g. IZheng et all 120051: iKravtsov erall I2004D 
show that the HOD is well approximated by a powerlaw 
Ns{M) oc M^ with e = 1, which is cut off below a certain 
minim al halo mass. Rath e r than this steep cut off, we 
follow IMandelbaum et aLl ()2005bi ) and assume a more 
gradual transition, and use e = 2 for halo masses smaller 
than Mchar, whilst e = 1 for halo masses larger than Mchar, 
where Mchar = 3M/j. Mh is the typical halo mass of a 
certain set of galaxies (for example the galaxies selected 
in a luminosity bin). The amplitude is determined by 
normalizing to the total number of satellites in the set. 
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4.1. Lensing signal from the halo model 

We now proceed to explain how the lensing signal is com- 
puted. The ensemble averaged tangential shear is the sum 
of the signal around central galaxies and satellites, since we 
cannot distinguish between them. We compute each contri- 
bution separately, starting with the signal around central 
galaxies. It is assumed that the central galaxies are located 
at the centre of the dark matter haloes. Two terms con- 
tribute to the lensing signal around central galaxies: the sig- 
nal coming from the halo where the galaxy resides (7t^cont)i 
and the signal from nearby haloes (7^ cent)- Hence the total 
signal around central galaxies is given by 



7i:Ccnt 



/i,cont ' 'i,ccnt' 



(15) 



The density profiles of the central haloes are assumed to 
be NFW, whi ch we compute usi ng the mass-concentration 
relation from iDuffv et al.l ()2008( ) given by Equation [71 By 
picking a central halo mass we can thus compute the tan- 
gential shear of the central halo term directly, as spectro- 
scopic redshifts are available for all lenses. 

The calculation of 7i^cent requires the power spectrum 
describing the correlation between the galaxy in the central 
halo and the dark matter of nearby haloes: 



P^::^,{k,Mn,T) = bg{Mh,T) 



PNL(fc) 

(27r)3 



Mii„ 



X / dvf{v)b{iy,r)ydmik,M), 

Jo 

with bg{Mfi,r) the bias of the ce ntral galaxy P^hik ) the 
non-linear power spectrum from ISmith et al.l (I2003D . and 
Udmik, M) the radial Fourier transform of the central halo 
density profile divided by mass: 



ydmik,M) 



M .,0 



dr47rrVd.n(r,A/)^^, (17) 
kr 



which we calculate using the analytical formula given in 
iPielorz et all (|2010l ). 

The dark matter profiles of adjacent haloes cannot over- 
lap, which is prevented by implementing halo exclusion. 
Different approaches to halo exclusion h a ve bee n used in the 
literature. For example, ICacciato et al.l ()2009f ) set the two- 
halo correlation function to zero below rigo , which leads to 
a sharp t runcation in t he ha lo models. We follow the ap- 
proach of [Tinker et al.l (|2005() , which leads to a more natu- 
ral smooth cut-off: the integral in Equation [12] is cut off for 
masses greater than Mum which is chosen such that the r2oo 
of the central halo does not overlap with the r2oo of nearby 
haloes: r2oo{Mh) + ?'2oo(Mim) = r. It should be noted that 
this choice, as any other halo exclusion approach, is an ap- 
proximation. Ultimately, numerical simulations should be 
used to provide improved estimates for P(.&nf 

The contribution of the satellites to the lensing signal 
consists of three terms: the signal from the subhalo where 
the satellite resides {'jl™^'^), the signal from the central halo 
in which the subhalo resides (7t\sat)i ^^'^ the signal from 
nearby haloes (7^ gat)- Hence the total signal around satel- 
lites is given by 

7t,sat = 7t,sat + 7t,sat + 7t,saf (1^^) 

First we compute the lensing signal of the subhalo, ^l'^^^'^, 
following iMandelbaum et aLJ (l2005bf ) . The density profile 



is assumed to follow an NFW profile in the inner regions. 
The outer regions of the subhalo arc tidally stripped of its 
dark matter by the central halo. Due to this stripping the 
lensing signal is proportional to r~^ at radii larger than the 
truncation radius. Ba sed on good agreem ent with numerical 
simulations. Mandelb aum et al.l (|2005bt ) chose a truncation 
radius of 0.4r2oo, and we use the same. This choice corre- 
sponds to roughly 50% of the dark matter being stripped 
from the subhalo. 

To compute the lensing signal induced by the halo where 
the subhalo resides, we calculate the power spectrum de- 
scribing the correlation between the subhalo and the dark 
matter profile of the central halo: 



P..tik,Mh 



1 



(27r)3fi 



diyf{iy)N,iM,Mh) 



(19) 



xydm(fc,M)yg(fc,M), 



with n the mean galaxy number density, which can be de- 
termined using n — p ^ dvf{v)—^^^-^ — —^ and j/g the ra- 
dial Fourier transform of the radial distribution of satellites 
around the central halo. We assume that the radial distri- 
bution of satellites follows an NFW profile with a concen- 
tration c^, g iven b y the mass-concentration relation from 
iDuffv et all (J2008D . To asses the sensitivity to the shape of 
the radial distribution of the satellites, we also calculate the 
(16) 7^sat term using a Cg that is varied by a factor of two. We 
find that this change mainly impacts the model signal at 
small scales: for a larger (smaller) concentration, the signal 
increases (decreases). At scales larger than a few hundred 
kpc, the change of the model signal is negligible. When we 
fit these adjusted models to the data, we find that the best 
fit model parameters do not change significantly. We con- 
clude that the signal-to- noise of our data currently does not 
enable us to discriminate between halo models with differ- 
ent radial distributions of satellite galaxies. 

Finally we compute the contribution from nearby haloes 
to the lensing signal around satellite galaxies: 



Pi'j(fc, M^, r) = %4? / '■" dvf{v)h{v, r)y^„,{k, M) 



(27r)3 / 

xf/d./(J)K.,0^^^^^^y.(fc,M^ 

(20) 

The three power spectra are converted into their respective 
shear signals using Equation [SJ and the contributions from 
the central galaxies and satellites are combined to yield 



7t = (1 - a) 7t,cont + a Jt,sa.t, 



(21) 



where a is the fraction of satellites of the sample. The 
resulting model is compared to the data. 

The lens sample is selected to cover a range in an ob- 
servable, such as luminosity or stellar mass, as the relation 
between the mean observable and the lensing mass is a use- 
ful constraint for simulations. The dark matter haloes of the 
lenses from such a sample have different masses, however, 
and it is therefore important to account for the scatter in 
the observable-halo mass relation. If the halo mass distri- 
bution is well-known, this can be done by integrating the 
models over the distribution of halo masses. Unfortunately, 
the distribution is generally not accurately known as the 
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lenses span a considerable range in observable, redshift and 
environment. A simpler approach is to study how the lens- 
ing mass is related to the mean halo mass for a given halo 
mass distribution. This a pproach, which was proposed by 
iMandelbaum et aP (J2006[) . provides the leading-order cor- 
rection for the scatter, and we use it in this paper. 



5. COMPARISON WITH DYNAMICAL MASS 

The dynamical mass traces the gravitational potential 
of a galaxy at small scales, and typically provides estimates 
of the total mass enclosed by the effective radius, which is 
of the order of a few kpcs. Comparison to the mass derived 
from strong lensing shows that both esti mates agree well for 
early- type galaxies (jBolton et al.l 120081 ). In contrast, weak 
lensing traces the gravitational potential at much larger 
scales, and the mass is usually determined within r2oo, 
whose values range between a few tens to a few hundreds 
of kpc. To study how the dynamical mass is related to the 
weak lensing mass, we measure the lensing signal for galax- 
ies divided into seven dynamical mass bins, as detailed in 
Table [1] The lensing signal of the stacked galaxies in each 
bin is shown in Figure \E[ We fit our halo model to the lens- 
ing signal in the distance interval between 50 kpc and 2 
Mpc. At scales smaller than 50 kpc the lensing signal is 
very noisy, since we do not have many sources at small sep- 
arations, and lens light contamination might bias the shear 
signal. At scales larger than 2 Mpc we measure the lens- 
ing signal using mainly sources that reside at the edge of 
the images, where the PSF ellipticity is large for the data 
taken prior to a change in the MegaCam configuratior|j(up 
to 15%), and the residual PSF systematics noticeably bias 
the lensing signal. We fit for the central halo mass and the 
satellite fraction, and use Eauation ll3l to compute the bias 
because the lensing signal is not well constrained at scales 
> 3Mpc. 

We impose two priors on the fits. Firstly, we do not fit 
halo masses that are lower than the mean stellar mass of 
the galaxies in the bin. This prior could introduce a bias 
if the assumed IMF is significantly different from the true 
one, leading to stellar mass estimates that are too high, 
but this is not expected to be the case. The second prior 
we impose is on the satellite fraction, which is not well con- 
strained by the data for the most massive galaxies and is 
anti-correlated with the best fit halo mass (see Appendix 
Ofor details). To prevent this from biasing the halo mass 
low, we limit the range of fitted satellite fractions to be less 
than 20% in the three highest dynamical mass bins as they 
contain galaxies that are expected to be nearly exclusively 
centrals. The best fit halo model for each bin is also shown 
in Figure IHl We find that the model fits the data well. The 
resulting best fit halo masses for the early- and late-type 
galaxies are shown in Figure [71 and detailed in Table [T] The 
error bars on the best fit halo mass (satellite fraction) in- 
dicate the la deviations determined by marginalizing over 
the satellite fraction (halo mass). 

For the early-type galaxies, we find that the dynam- 



^ In November 2004, the lens L3 was accidentally mounted 
incorrectly after the wide-field corrector had been disassembled. 
As this surprisingly led to a significant improvement in the image 
quality for the u*-, gi'-,and r'-band, the new configuration was 
kept. About 20% of the RCS2 survey was obtained prior to this 
change. 
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Fig. 7. The best fit halo mass as a function of the mean 
dynamical mass. The red squares (blue triangles) denote 
the halo mass for the early- types (late- types). The early- 
/late-type division is based on the brightness profiles of the 
lenses. The dynamical mass correlates well with the lensing 
mass for the early-type galaxies, but not for the late-type 
galaxies. 



ical mass correlates well with the halo mass. The halo 
mass is ^10 times larger than the mean dynamical mass 
for Mdyn < 1 X 10^"'^ M0, which increases to a factor ~50 
for the highest dynamical mass bins, as the galaxy dark 
matter haloes extend far beyond the effective radius. To 
establish whether we can scale the dynamical mass to the 
lensing mass, we replace Re with the best fit lensing r2oo 
in Equation [TJ We find that the rescaled dynamical masses 
are 8 times larger than the best fit lensing masses for Dl 
and D2, but the difference decreases for the more massive 
bins: the rescaled dynamical mass is only 40% larger than 
the best fit lensing mass for D7. We therefore cannot simply 
rescale the mean dynamical mass to the lensing mass. Note 
that at the high mass end, galaxies predominantly live in 
groups and clusters. With lensing we fit the halo mass of 
the entire structure, whereas the dynamical mass is deter- 
mined for the individual galaxy only. 

We observe that for the late-type galaxies the halo mass 
does not correlate well with the mean dynamical mass. In 
particular, the best fit halo masses of the D5 and D6 late- 
type bins are low. These low values may be explained if 
rotation constitutes a major part of the observed velocity 
dispersions of late-type galaxies, leading to an overestima- 
tion of the dynamical mass. Additionally, the effective ra- 
dius for some late-type galaxies at high redshift may be 
overestimated, since a significant fraction consists of multi- 
ple objects with small separations as we observed in Section 

m 

For early-type galaxies the dynamical mass is a useful 
tracer of the total mass at small scales, but it appears to 
be less reliable for late-type galaxies. How the dynamical 
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Fig. 6. The lensing signal AS for each dynamical mass bin as a function of physical distance from the lens. The green 
dashed line shows the Jt'c^nt term, the blue dashed line the 7t^ccnt term, the green dotted line the 74™'^ term, the red 
dotted line the 7(^sa,t term, the blue dotted line the 7t^sat term, and the black line shows the sum of the terms. The 7^^ g^j 
term causes a prominent bump for the two lowest dynamical mass bins, which indicates that a significant number of 
lenses in these bins are satellites. 

Table 1. The dynamical mass results 



Sample 


log(Mdyn) 


?^lcns 


(^) 


(Mdyn) 


/late 


^^cat-ly 


^carly 


Mt'" 


^lato 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Dl 


[10.00, 10.50] 


2011 


0.08 


1.96 


0.44 


2-00+lil 


24+°" 

"■^^-0.09 


4.9Q+tit 


0-i7+Joll 


D2 


[10.50, 11.00] 


4 752 


0.10 


5.91 


0.35 


3.98tlil 


r, .1+0.07 

U.4i_Qog 


0-74tlf3 


n 11+012 


D3 


[11.00,11.25] 


2 762 


0.13 


13.2 


0.25 


17.8tr.50 


r, 14+0.07 


4.27l|:?| 


0.191°:!^ 


D4 


[11.25,11.50] 


2 281 


0.16 


23.6 


0.16 


28--2tlil 


n qi+009 
U.Oi_o.08 


31.6t}^i 


0.001°:°^ 


D5 


[11.50,11.75] 


1715 


0.22 


41.7 


0.07 


W5tll:l 


o.2o«:°° 


8.91_ggo 


o.ool°:gg 


D6 


[11.75, 12.00] 


935 


0.32 


72.5 


0.05 


295+?^ 


o.2ot°:°g 


3A7tt:,l 


o.oot°:^° 


D7 


[12.00, 12.50] 


380 


0.39 


137.4 


0.07 


^GSlilll 


o.2o«:°° 


219tllt 


0.201°:°° 



Notes. (1) the dynamical mass range of the bin; (2) the number of lenses; (3) the mean redshift; (4) the mean dynamical mass in 
units of 1O^°M0; (5) the fraction of late-type galaxies; (6) the best fit halo mass for the early-types in units of 10^^h~^MQ; (7) 
the best fit satellite fraction for the early-types; (8) the best fit halo mass for the late-types in units of 10^^ h~^ Mq; (9) the best 
fit satellite fraction for the late-types. 



mass changes for galaxies where rotation is important, or 
for galaxies that are populated over a large range of red- 
shifts, may be studied with numerical simulations. In any 
case, it is not clear how to translate a dynamical mass esti- 
mate into a total mass estimate of the halo of a galaxy. With 
weak lensing we measure the total halo masses of galaxies 
directly, providing estimates that can easily be compared 
to simulations. 



6. LUMINOSITY RESULTS 

The optical luminosity is a readily measured quantity 
which is related to the stellar mass, and hence the baryonic 
content of a galaxy. Therefore, we continue by measuring 
the lensing signal as a function of luminosity. We divide our 



lens sample into eight luminosity bins, as detailed in Table 
[5] We measure AE of the stacked lenses and show the 
results in Figure HI together with the best fit halo model. 
The amplitude of the lensing signal clearly increases for 
the brighter galaxies as expected. Furthermore, the shear 
from the 7j^sat term causes a prominent bump for the 
fainter lenses, but not for the brighter ones. This indicates 
that a considerable fraction of the low luminosity lenses 
are satellites. We split the lenses into early- and late-types 
using the frac-deV parameter as before, and study the 
signals separately. 



There are various issues we have to address before we 
can interpret the measurements. First of all, lens galaxies 
scatter between luminosity bins due to luminosity errors. 
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Fig. 8. The lensing signal AS for each luminosity bin as a function of physical distance from the lens. The green dashed 
line shows the 7(^cont term, the blue dashed line the 7t^cont term, the green dotted line the 74^3™'^ term, the red dotted line 
the 7t^sat term, the blue dotted line the 7t^sat term, and the black line shows the sum of the terms. A significant fraction 
of the low luminosity lenses are satellites in larger haloes, as the jl'lat term causes a prominent bump at ~ IMpc in the 
lensing signal. 

Table 2. The luminosity results 



Sample 


Mr 




"Ions 


(^> 


/late 


(^learly^ 


^,^early 


yearly 


T early 
-^200 


{L[-'^) 


M]f'' 


late 

a 


rlate 
-1^200 




(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


LI 


-21.0,- 


-20.01 


3 563 


0.08 


0.44 


1.17 


4.78t^;S 


0.37«;S? 


1.32±0.01 


1.23 


1.66±?-I 


Q.m+_l-l\ 


1.26±0.01 


L2 


-21.5,- 


-21.0] 


2 772 


0.10 


0.35 


2.14 


3.67+lil 


0.38t°:g«8 


2.36±0.01 


2.34 


3.601*:^^ 


o.o6l^:jo 


2.52±0.01 


L3 


-22.0, - 


-21.5] 


3 064 


0.13 


0.29 


3.24 


16.5tlil 


o.i9t°:°? 


4.01±0.02 


3.63 


2.14l«^ 


0-35_o;i2 


3.67±0.01 


L4 


-22.5,- 


-22.0] 


2 370 


0.16 


0.20 


5.01 


91 q-t-7.58 
^1.0„7^28 


O.S4t°°ol 


6.57±0.03 


5.62 


16.6±^°5? 


o.ootH^ 


6.50±0.04 


L5 


-23.0, - 


-22.5] 


1658 


0.20 


0.11 


7.59 


105t^i-9 


0-05trol 


12.3±0.1 


8.91 


12.7+1*-? 


o.2ol«:?o 


lO.OiO.l 


L6 


-23.5,- 


-23.0] 


1453 


0.32 


0.03 


11.0 


267t«i-i 


o.i9«:?^ 


21.2±0.2 


13.8 


14lti« 


0.20j:«:~ 


23.4±1.7 


L7 


-24.0, - 


-23.5] 


607 


0.41 


0.05 


15.8 


570tI?o' 


o.oo«:g« 


45.2±0.8 


21.9 


306^^7^ 


o.oot«:go 


- 


L8 


-24.5,- 


-24.0] 


83 


0.47 


0.04 


22.9 


818lS^ 


O-OOt'ofo 


68.0±3.6 


- 


- 


- 


- 



Notes. (1) the magnitude range of the bin; (2) the number of lenses; (3) the mean redshift; (4) the fraction of late-type galaxies; 
(5) the mean luminosity for the early-types in units of 10^'^ Lq\ (6) the best fit halo mass for the early-types in units of 10^^ h~^ Mq\ 
(7) the best fit satellite fraction for the early-types; (8) the total luminosity within r2oo for the early-types in units of 10^^ Lq; (9) 
the mean luminosity for the late-types in units of lO^^Lo; (10) the best fit halo mass for the late-types in units of \Q^^h~^MQ\ 
(11) the best fit satellite fraction for the late-types; (12) the total luminosity within r2oo for the late-types in units of lO^^L© . 



If the luminosity errors are large compared to the width 
of the bins this could potentially introduce a bias. This 
bias is greatest at the highest luminosities, where the 
luminosity function is steep. In this case, on average more 
low luminosity (and mass) galaxies scatter into the higher 
luminosity bins, biasing the best fit halo mass low. The 
average absolute magnitude error is ^-^0.03 for z < 0.33, 
and ~0.07 for z > 0.33, small compared to the minimal 
bin- width of 0.5. We find that the induced bias is relevant 
for the L7 and L8 bins of the early-types only, with 
corrections of 4% and 7% respectively. The corrections are 
smaller than the measurement errors on the halo mass 
for these bins. We detail the calculation of the correction 



factor in Appendix \K\ 

When we fit a halo mass to the stacked shear signal 
of galaxies within a luminosity bin, the resulting mass 
is not equal to the mean halo mass, nor to the cen- 
tral m ass of the origina l distrib u tion (iTasitsiomi et al.l 



2004 iMandelbaum et all l2005bl: ICacciato et al.l 120091: 
Leauthaud et al.l l2010l ) because the distribution in halo 



mass is not uniform (in addition, the NFW profile itself 
depends on mass). It is useful to convert the measured 
lensing mass to the mean halo mass to allow comparison 
with simulations. The correction we have to apply depends 
not only on the distribution of halo masses for a given 
luminosity, but also on the halo mass function. Since the 



12 



Edo van Uitert et al.: The relation between baryons and dark matter 



^ 10'' 

k 10'^ 

o 
o 

CN 

^ 10" 

600 

" 500 

^ 400 

^^ 300 

g 200 

100 




(q) 



1-1 



*"'M--'] 



I 



(b) 



J .. 



i 



■> ^ n 



10' 



L. [Lo] 



10' 



Fig. 9. The best fit halo mass (top), and the mass-to-light 
ratio (bottom) as a function of mean luminosity. The red 
squares (blue triangles) denote the early-type (late-type) 
results. The division in early-/late-types is based on the 
brightness profiles of the lenses. The dashed lines are the 
power law fits, with values as indicated in the text. 



halo mass function is a declining function — steeply at the 
high mass end — we will preferentially select lower mass 
haloes. Hence, the underlying function from which we draw 
our galaxies is the halo mass function convolved with the 
halo mass distribution. In Appendix [B] we discuss how we 
calculate the correction factor that we apply to obtain the 
mean of the halo mass in each luminosity bin. The values 
are given in Table [BTI and range between 5-30%. 

The best fit halo mass for each luminosity bin, cor- 
rected for the scatter and the width of the halo mass 
distribution, is given in Tabled and is shown as a function 
of luminosity in Figure IH^. The error bars on the halo 
masses are the Icr deviations determined by marginalizing 
over the satellite fraction. We fit a powerlaw of the form 
M200 = Mo,l{L/LoY^, with a pivot io = lO"Lr,0. As 
the errors of the best fit halo masses are asymmetric 
due to the constraints we impose on the halo model fits, 
we fit the powerlaw directly to the shear measurements 
(with symmetric error bars). Hence we do not fit for 
the halo mass for each bin, but determine the best fit 
Mo,L and Pl for all bins simultaneously, whilst fitting the 
satellite fraction for each bin separately. Note that the 
best fit satellite fractions from this approach are close 
to the values given in Table [31 For the early-types, we 
find Mo,L = l-93i°:i| x lO^^/i-^M© and Pl = 2.341°°^, 
and for the late-types Mq^l = 0.43J:o:i? x IO^^H-^Mq 



and Pl = 2.2^gg, as shown in Figure ^j^. The error on 
Mo,L [Pl) is determined by marginalizing over Pl (Mq^l). 
We' show the 67.8%, 95.4% and 99.7% confidence hmits 
of the two powerlaw fits in Figure [TUl The results for 
the early-types are better constrained because we have 
more early-type galaxies in our lensing sample. These are 
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Fig. 10. The 67.8%, 95.4% and 99.7% confidence limits of 
the powerlaw fits between luminosity and halo mass, in 
red (blue) for the early-type (late-type) galaxies. The red 
square (blue triangle) indicates the best fit of the early- 
types (late-types). The powerlaw fit for the early- types is 
better constrained than for the late-types, because the ma- 
jority of galaxies in our lens sample are early- types. The 
early- types also reside in more massive haloes, and conse- 
quently produce a stronger lensing signal. 



also more massive than the late-type galaxies and hence 
produce a stronger lensing signal. 

We compare our an a lysis t o two previous weak lensing 
studies. iHoekstra et al.l (J2005t) measured the lensing signal 
of ~ 1.4 X 10^ isolated galaxies with photometric redshift 
0.2 < z < 0.4 in the RCS. In the i?-band, they found 
a virial mass of M^°^ = 7.5til x IO^/i-^Mq for a 
galaxy of luminosity Lfj = IQ^^h'^'^LQ, and a powerlaw 
index of P™^ ~ 1.6 ± 0.2. We use the transformations 
from Lupton (2005)0, and find that r k, R + 0.24 for 
the early-type galaxies in our sample, which make up 
the majority of the lenses. We convert Ln to L^, use 
our powerlaw fit to predict A'/2oo, and convert that to 
the virial mass by increasing it by 30%. We find that 
Afvir = (7.2 ± 1.5) X 10"/i-iA/o for a Lb. ^ 10| " /i^^Lff 
galaxy, in good agree ment with iHoekstra et al.l (|2005[ ). 
The powerlaw index of IHoekstra et all (120051 ) is shallower 
than the Pl = 2.34+Q;5g that we find. A possible expla- 
na tion is that a f ractio n of the low luminosity galaxies 
in IHoekstra et al.l (|2005f ) are satellites, whose masses are 
biased high due to the added lensing signal of nearby 
galaxies, flattening the powerla w index . We note two 
caveats: the lens sample of Hockst ra et al.l (J2005l ) does not 
exclusively consist of early-types, and the lens samples we 



* http://www.sdss.org/dr7/algorithms/sdssUBVRITrans 
form.html#Lupton2005 
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co mpare reside in diffe r ent en vironments. 

iMandelbaum et all (|2006[) present results for 3.5 x 10^ 
galaxies using SDSS data. Galaxies arc divided into 
early-types and late-types based on their brightness profile 
(using the same selection criterium that we have applied to 
our lenses), and are studied in bins of absolute r-band mag- 
nitude. To compare the results, we conv ert our luminosities 
accord ing to the definitions used in IMandelbaum et al.l 
(|2006f ): the absolute magnitude is calculated using a 
k-correction to z=0.1, the distance modulus is calculated 
using h = l.O and a passive evolution term is included 
which is given by 1.6(z — O.f). As a result, we decrease 
the absolute magnitudes of our lenses by roughly one 
mag nitude. Additionally , we in crease our masses by 30% 
since IMandelbaum et aLl ()2006l) define the halo mass using 
180p instead of 200pc- There are various other differences 
between the analyses, such as the use of a different correc- 
tion factor for the width of the halo mass distribution, a 
different cosmology, a different mass-concentration relation 
for the NFW profiles, and differences in the modelling of 
the Icnsing signal. These differences are expected to have a 
minor impact on the best fit halo mass, but they limit the 
accuracy of a detailed comparison. 

Matching our luminosity bins to those of 
IMandelbaum et ah! (|2006[) closest in mean luminosity, 
we find that the best fit halo masses for the early- and 
late-type galaxies are generally in agreement. To quantify 
whether the results are consistent, we fit a powerlaw of the 
form Afiso = Mo,L(i/io)'^^, where Lq = 1.2 x IO^^H-'^Lq. 
The tilde indicates that th e lumi nosity is calculated 
following IMandelbaum et aP (J2006f ). The powerlaw is 
fitted to the best fit halo mass directly, and the weights 
of the measurements are calculated from the error bars 
through which the model passes, i.e., if the model is 
larger (smaller) than the data point, we use the pos- 
itive (negative) error bar. For the early-types we find 
Mo,L = 7.3tl;\ X 1 0"/i-iM,T, and 0f = 2.1 ± 0.2 for our 
data, while using IMandelbaum et all ()2006D results we 
find Mo,L = ll-StiJ X lO"/i-iM0 and ^i = 2.3 ± 0.2, 
in fair agreement with our findings. For the late- types 
we find Mo.L = 2.7t?;8 ^ lO^/f^Mp, and I3r = 3.0 +^g, 
while using the results of IMandelbaum et aLl (|2006t ) we 



find Afo L = 7.8 ± 1.1 X W^h-^Mt^ and 



1.1 



+0.3 
-0.4- 



The 



results from ^M andelbaum et all (|2006l ) prefer a shallower 
slope and a higher offset, but the fits are consistent. 



6.1. Mass-to-light ratio 

A large number of the galaxies in our brightest luminos- 
ity bins reside in groups or small clusters. To identify 
those lenses, we cross-correlate our lens sample with the 
preliminary RCS2 cluster catalogue, to be presented in 
a future publication. We take galaxies with a separation 
< 250ft-~^kpc from the cluster centre, and within 0.05 from 
the cluster redshift, to be cluster members. Using these cri- 
teria, we find that from L5 to L7, 3%, 26%, 43% of the 
late-type galaxies, and from L5 to LB, 12%, 31%, 48% and 
66% of the early-type galaxies can be associated with clus- 
ters. The best fit halo mass of these galaxies is the mass 
of the group or cluster within r2oo, while the luminosity 
is only measured for the lens galaxy. The resulting mass- 
to-light ratio, shown in Figure [Hb, is therefore higher than 



what we would measure for the individual galaxies, or for 
the clusters. 

To obtain the mass-to-light ratios of the groups and 
clusters, wc estimate the amount of additional luminosity 
coming from other cluster members within 7'2oo- We as- 
sume that the spectral energy distributions (SEDs) of the 
galaxies physically associated with the lens are similar to 
the SED of the lens, and convert their apparent magni- 
tudes to absolute magnitudes using the same conversion 
that has been used for the lenses. The apparent magni- 
tude s we use are those f rom the photometric catalogues 
from lGilbank et al]()2011h . As these catalogues do not cover 
all fields (e.g. the fields in the uncompleted patch 1303), 
only ~90% of the lenses are used for the calculation of ^200- 
We measure the source galaxy overdensity as in Section [231 
using all the galaxies with ?niow < mr < 24, where toiow is 
the magnitude of the brightest galaxy that resides at the 
lens redshift, and calculate the mean luminosity overdensity 
as a function of lens-source separation, miow is determined 
by selecting the brightest galaxy in th e photometric red- 
shift catalogues from Ilbert et al.l ()2006l ) that resides at the 
redshift of the lens or higher. We sum the luminosity over- 
density to r2oo and add it to the lens luminosity to obtain 
the total luminosity within r2oo, ^200- To make sure that 
we do not miss a signicant fraction of L200 from galaxies 
with rrir > 24, we also calculate L200 using an upper limit 
of 23.5, and find that the results do not change significantly. 
The values of L200 are given in Table j^l We show the mass- 
to-light ratio A/200/^200 a-s a function of £200 in Figure [TTJ 
For ^200 < lO^^L© we calculate the weighted mean, and 
find a value of M2oo/-£20o = 42± 10 for early-type galaxies, 
whilst M2oo/i2oo = 17± 9 for late-type galaxies. The total 
mass-to-light ratio increases with L200 for the early-types 
to ^--^180 at L200 = 5 X IO^^Lq. The total mass-to-light ra- 
tio is roughly a factor of two larger for early-types than for 
late-types. This suggests that the difference in the best fit 
halo mass between early- and late-types for a given lumi- 
nosity is not solely due to the fact that early-types reside 
in denser environments, but is at least partly intrinsic. The 
value of L200 for the L7 late-type bin could not be robustly 
determined, and is excluded from the results. 

We com p are ou r results to the A/200 /-^200 from 

ISheldon et al.l (|2009al ibl) which have been determi ned for 
the clusters in the maxBCG catalogue ( Koester et al.l 
I2OO7D . The quoted values of L200 in their work have been 
measured in the i-band, and are calculated using a k- 
correction to z = 0.25. We convert them to the r-band 
luminosities we use by accounting for the mean difference 
between i-band and r-band absolute magnitudes of early- 
type galaxies at z = 0.25, the mean difference between the 
k-corrections to z = 0.25 and z = 0.0, and the difference 
between the z-band and r-band solar magnitudes. The final 
conversion factor is small as the corrections partly cancel 
each other, and we convert their luminosities to our defi- 
nition by multiplying them by 1.06. Note that we do not 
account for differences in the redsh ift evolution of the lu- 
minosities, as it is not mentioned in ISheldon et al.l ()2009a[) 
which correctio n, if any, they have u sed. The converted 
^•^20o/-^20o from ISheldon et ahl ([20093) are indicated with 
the hatched area in Figure[TTJ The mass-to-light ratios over- 
lap, and the ratios we have determined, for individual galax- 
ies at low luminosities, and for galaxy groups and small 
clusters at high luminosities, are naturally extended to the 
^£200/^^200 of clusters from the maxBCG cluster sample. 
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Fig. 11. The mass-to- light ratio using the total halo mass 
and luminosity within ^200, as a function of £200- The 
red squares (blue triangles) denote the early-type (late- 
type) results. The hatched area indicates the converted 
M200/L 200 of the maxBCG clusters from ISheldon et al.l 
(|2009af ). The M2oo/-^20o for individual galaxies at low lu- 
minosities are naturally extended to the ratios for the 
maxBCG clusters. 



6.2. Satellite fraction 

Figure [T^] shows the best fit satellite fraction as a function 
of luminosity. The satellite fraction is decreasing with in- 
creasing luminosity for the early-type galaxies, from ~ 40% 
at Lr = 10^°Lq to < 10% at L^ = lO^L©. For the late- 
type galaxies, no clear trend with luminosity is observed, 
and the satellite fraction has a value of 0-20%. The satel- 
lite fractions are not well constrained for the highest lu- 
minosity bins. As demonstrated in Appendix O the sum 
of the halo model satellite terms has the same shape as 
the central term at the high halo mass end. As a result, 
the halo model fit cannot discriminate between the two 
profiles. The implementation of a more sophisticated de- 
scription of the truncation of the subhaloes is necessary 
to improve the constraints on the satellite fraction at the 
hi gh luminosity/st e llar nr ass end. For instance, recent work 
bv lLimousin et al] ()2009() suggests that massive early-type 
satellite galaxies are stripped of a far larger fraction of their 
dark matter than the 50% we have assumed so far, and we 
dis cuss the implications in A ppendix [Cl 

iMandelbaum et aP (|2006f ) find a satellite fraction of 10- 
15% for late-type galaxies, independent of stellar mass or 
luminosity. The satellite fraction for early-types decreases 
with luminosity from 27% at (L/Lq) = 1.1 to 15% at 
{L/Lq) = 4.9, and both trends are consistent with our find- 
ings. 



Fig. 12. The best fit satellite fraction as a function of 
mean luminosity. The red squares (blue triangles) denote 
the satellite fraction for the early-type (late-type) galax- 
ies. The satellite fraction decreases with luminosity for the 
early-types, and no trend is observed for the late- types. The 
dashed area indicates the area excluded by the prior on the 
satellite fraction. 



7. STELLAR MASS RESULTS 

The stellar mass of a galaxy is believed to be a better tracer 
of the baryonic content of a galaxy than the luminosity, 
as it is less sensitive to recent star formation. Therefore, 
we divide our lens sample into seven stellar mass bins and 
study the lensing signal. The details of the samples are 
listed in Table [31 Figure [T3] shows the lensing signal of the 
stacked lenses in each bin, together with the best fit halo 
model. Similar to the luminosity results, we find that the 
lensing signal increases with stellar mass, and observe the 
presence of the ^t'lat bump for the lower stellar mass bins. 
We split the lens sample into early- and late-types using 
the frac-deV parameter as before, and study the signals 
separately. 

To interpret the results, we have to account for a 
number of issues. The random stellar mass errors are 
about 0.1 dex, independent of stellar mass, and do not 
include the systematic error. The random error determines 
the scattering of lenses amongst bins, and its value is 
large compared to the bin width. We calculate the bias 
resulting from this scatter in Appendix |21 and find that 
the best fit halo masses have to be corrected with a factor 
ranging between 0.9 — 1.4. Once corrected for the scatter, 
we convert the lensing mass to the mean halo mass. This 
correction has already been introduced in Section [51 and 
we discuss in Appendix [BJ how wc calculate it. We increase 
the corrected halo mass accordingly to obtain the mean 
halo mass (see Table [BTJ for details). 

The resulting halo masses are given in Table [21 and 
shown in Figure 1141 This figure shows that the relation 
is different for early- types and late- types. Below a stellar 
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Fig. 13. The lensing signal AS for each stellar mass bin as a function of physical distance from the lens. The green 
dashed line shows the 7t^cont term, the blue dashed line the 7^cont term, the green dotted line the 7t'sat'^ term, the red 
dotted line the 7(^sat term, the blue dotted line the 7^sat term, and the black line shows the sum of the terms. Similar to 
the lensing signal of the luminosity bins, we find that the jl^at term causes a clearly noticeable bump at ^ IMpc in the 
lensing signal for the low stellar mass bins, which indicates that a significant fraction of these galaxies are satellites. 

Table 3. The stellar mass results 



Sample 



log(M.) 
(1) 



^Icns 

(2) 



(3) 



(M.> /late 



(6) 



early 



(7) 



a /foarly 
^'•-'*,200 

(8) 



(9) 



(10) 



(11) 



SI 


[10.00, 


10.50] 


3 359 


0.08 


2.04 


0.51 


1 qr+2.30 


o.44l°:g^ 


S2 


[10.50, 


11.00] 


5 870 


0.11 


6.03 


0.28 


8.39+2I1 


0.191H1 


S3 


[11.00, 


11.25] 


2 428 


0.15 


13.2 


0.10 


i4n+6-2i 


n 44+009 

U.44_g Qg 


S4 


[11.25, 


11.50] 


1631 


0.20 


24.0 


0.05 


1351-- 


o.oolH^ 


S5 


[11.50, 


11.75] 


1505 


0.34 


41.7 


0.03 


400l^°^e 


0.20t»;~ 


S6 


[11.75, 


12.00] 


396 


0.41 


69.2 


0.05 


640±^i? 


0.06±«;J^ 


S7 


[12.00, 


12.50] 


48 


0.48 


123 


0.02 


7221^31 


0.201°;* 



2.15±0.004 


0.56tSJ^ 


10+0" 


2.04±0.001 


7.66±0.02 


10 7+-'5-ii 


0.121^8 


6.75±0.01 


16.8±0.1 


1.49li^i 


o.4olg:« 


14.4±0.1 


40.3±0.3 


2M+l'i 


0.701°:^,? 


27.3±0.4 


88.9±1.1 


65.21^^2 


0.20t°;°° 


71.0±4.6 


230±6 


001 +353 

ZZl_22i 


0.20lg;~ 


- 


385±28 


- 


- 


- 



Notes. (1) the stellar mass range of the bin; (2) the number of lenses; (3) the mean redshift; (4) the mean stellar mass in units of 
10^'^ Mq; (5) the fraction of late-type galaxies; (6) the best fit halo mass for the early-types in units of IO'^^/i'^Mq; (7) the best 
fit satellite fraction for the early-types; (8) the total stellar mass within r2oo for the early-types in units of IO'^'^Mq; (9) the best 
fit halo mass for the late-types in units of 10^^h~^MQ; (10) the best fit satellite fraction for the late-types; (11) the total stellar 
mass within r2oo for the late- types in units of 10^° M© . 



mass of IO^^Mq, the halo mass is similar for both galaxy 
types, but for stellar masses larger than 10^^ Mq the halo 
masses of early-type galaxies are more massive for a given 
stellar mass than the halo masses of late-type galaxies, 
and increase more steeply with stellar mass. These trends 
in the stellar mass t o halo mass relation are i n agreement 
with those found by iMandelbaum et al.l (I2006D . 

We fit a powerlaw of the form M200 — 
Mo^m{M^/Mo)^", with Mq = 2 X 10"Mq, fit- 
ting the lensing measurements simultaneously as we 
did for the luminosities. For the early-types, we find 
Mo,M = 8.1 ± 0.6 X lO^^h-^MQ and Pm = 1-9 ± 0.1, 
and for the late-types Mq^m = '^-^tl'^ x IO^^/j-^Mq and 



(3m = 1.2 ± 0.4. These fits are shown in Figure [T3] as the 
dashed red and blue lines for the early- and late-type 
galaxies respectively. We show the 67.8%, 95.4% and 99.7% 
confidence limits of the two powerlaw fits in Figure [TSl 

In order to com pare with the results of 
IMandelbaum et ahl (|2006[ ). we lower their halo masses 
by 30% to account for the difference between Mvir and 
Af2oo- We compare the halo masses of the bins with 
comparable mean stellar mass, and find that the best 
fit halo masses generally agree well. We fit a powerlaw 
between stellar mass and halo mass to their results, and 
the dashed contours in Figure [TS] show the resulting best 
fit normalisation and slope. The results of the late-types 
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Fig. 14. The best fit halo mass as a function of mean stellar 
mass. The red squares (blue triangles) denote the early- 
type (late-type) galaxies. The separation of the lenses into 
early/late-types is based on their brightness profiles. The 
dashed lines are the power law fits, with values as indicated 
in the text. For stellar masses lower than IQ^^Mq the best 
fit halo masses of early- and late- type galaxies are similar, 
but for M* > IQ^^Mq we find that the best fit halo masses 
of early-types are greater. 



5a. ^ 




5x10 
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Fig. 15. The 67.8%, 95.4% and 99.7% confidence limits of 
the powerlaw fits between stellar mass and halo mass, in 
red (blue) for the early-type (late-type) galaxies. The solid 
contour lines result from fitting the powerlaw to all the lens- 
ing data as described in the text. The dashed contours are 
the results from fitting the powerlaw b etween stellar mass 
and h alo mass to the measurements in iMandelbaum et al.l 
(l2006l) . 



agree, although the errors are large. For the early-types, 
IMandelbaum et all ([2006) find a somewhat steeper slope 
and a higher offset. Note that this difference is mostly 
driven by their highest stellar mass bin, for which they fit 
a halo mass that is 50% larger than what we find for our 
corresponding bin. If we exclude that point from the fit, 
th e la contours o verlap . 

iMoster et al.l ()2010l ) used numerical simulations to 
predict the relation between stellar mass and halo mass. 
We find that for M, < 4 x 10^^ M©, the halo masses we 
have determined are about 1-2 a lower than their models. 
At higher stellar masses, the discrepancy is significantly 
larger. Not only the ir model, but also t he models of various 
other groups (e . g. Wang et al. 2006t Croton et al. 2006t 



Somerville et al.l 120081: iBehroozi et all 120101: iNeisteinet ahl 
20111) predict that the halo masses of galaxies with a stellar 



mass > 10 Mq increases rapidly as a function of stellar 
mass, a trend we do not observe in our measurements. This 
would imply that the predicted relation between stellar 
mass and halo mass for galaxies with A/* > 4 x 10^ ^Mq 
is too steep, possibly because the relation has not yet 
been well constrained by observations in this mass range. 
Although contamination of the high stellar mass bins 
by unresolved mergers may bias the best fit halo masses 
low, we estimate that this is not sufficient to explain the 
discrepancy. 



7.1. Baryon conversion efficiency 

To study the efficiency of star formation as a function of 
stellar mass, we measure the baryon conversion efficiency 
1] — Mf,/{Mh X /(,), where /(, = Q,b/Q,M is the cosmological 
baryon fraction. We cannot simply use the mean stellar 
and halo mass, because we measure the halo mass of the 
environment where the galaxy resides. The mean stellar 
mass, however, is determined using the individual lenses 
only, which leads to an underestimation of r]. To account 
for this, we estimate the additional amount of stellar 
mass within r2oo assuming that the SEDs of the cluster 
members are similar to that of the lens galaxy. Under that 
assumption we determine M<,_200 = (-M*) x (-^20o/-^r), 
where L200 is the total luminosity within 7^200 as discussed 
in Section 16.11 The error bars assume that the number 
of source galaxies in each radial bin follows a Poisson 
distribution. We give the values of M,^200 in Table [21 and 
plot 77 as a function of M*_200 in Figure [T6l 

The stars that make up the diffuse intracluster light 
(ICL) also contribute to the total stellar mass. The ICL 
typically makes up 10- 20% of the stellar light in galaxy 
groups and clusters (see iGiodini et al.ll2009l and references 
therein). We do not account for the additional stellar 
mass from the ICL, because our lens sample consists of 
a mixture of isolated galaxies and galaxies in groups and 
clusters. The average contribution from the ICL is hard to 
determine, particularly because the contribution for low 
mass structures is very uncertain. The ICL is expected 
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Fig. 16. The baryon conversion efficiency ry as a function 
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squares (blue triangles) denote the early-type (late-type) 
galaxies, rj is smaller for the early-types than for the late- 
types for M,,200 > 10"Mq. 



to be of importance for the S6 and S7 bins only, as they 
contain the largest fraction of cluster associated galaxies, 
and the derived values of rj might at most increase with 
10-20%. 

We find that rj decreases from ^40% at 
M* -- 5 X 10^°AfQ to a minimum of --10% for a 
stellar mass AfH.,200 = lO^^M©, and seems to increase 
again at higher stellar masses. The baryon conversion 
efficiency for the late- types is higher, and no clear trend is 
observable because of the large errors. For some bins r] is 
larger than unity, but the error bars cover the reasonable 
range of 77 < 1. The value of Mh._2oo for the S6 late-type bin 
could not be robustly determined, and is excluded from 

th e results. 

iHoekstra et al.l (J2005t) divide the lens sample in red and 
blue galaxies based on their B — V colour, and find that 
the baryon conversion efficiency for isolated blue galaxies 
in the magnitude range 18 < Re < 24 is about tvifice the 
value found for isolated red galaxies. Although we cannot 
compare the results in detail due to differences in the type 
selection and differences in the adopted IMF, our results 
also suggest a larger value for t] for late-type galaxies in 
the range M y. 200 > 10^^M(t,. A s i milar trend has also been 
observed in iMandelbaum et ahl (|2006[ ) for stellar masses 
M* > 10^^, but note that the baryon conversion efficiencies 
were determined using M* instead of M*,200, and the 
values are therefore lower limits. 



7.2. Satellite fraction 

In Figure [T7] we show the satellite fraction for the early- 
and late- type galaxies as a function of stellar mass. The 



Fig. 17. The best fit satellite fraction as a function of the 
mean stellar mass. The red squares (blue triangles) denote 
the early-type (late-type) results. The satellite fraction de- 
creases with stellar mass for the early- types, and no trend 
is observed for the late- types. The dashed area indicates 
the area excluded by the prior on the satellite fraction. 



satellite fraction of the late-types is only well determined 
for the SI and S2 bins, and appears to be constant as a 
function of stellar mass, with a value of ^10%. The satel- 
lite fraction of the early-types is 45% for the lowest stellar 
mass bin, but decr e ases t o < 10% for M, > 2 x IO^^Mq. 
IMandelbaum et aD (l2006t ) find a satellite fraction of about 
10-15% for late- type galaxies, in dependent of stellar mass or 
luminosity. For the early-types, IMandelbaum et all ()2006() 
find that the satellite fraction decreases with stellar mass 
from 50% at 10^° Mq to roughly 10% at 3 x IG^Mq, con- 
sistent with our findings. 

7.3. Dependence on redshift 

The stellar mass of a galaxy and the dark matter content 
of its halo evolve with time. The stellar mass increases 
as galaxies form stars and merge with satellites and other 
galaxies. Satellite galaxies residing in subhaloes are tidally 
stripped of their dark matter, whilst the dark matter con- 
tent of central haloes increases due to mergers. The evolu- 
tion of the relation between stellar mass and dark matter 
content of g alaxies has been studied with numer ical simu- 
lations fe.g. iMoster et al.ll2010l: IConrov fc Wechs ler 200^. 
These simulations predict that the dark matter content of 
haloes that host galaxies of M* > IO^^Mq increases faster 
than the stellar mass, while the stellar mass grows faster 
for haloes hosting galaxies of M* < 10^^ Mq. 

To study this, we bin the early-type galaxies in stellar 
mass and redshift, and measure their halo mass. To avoid 
the degeneracy between halo mass and satellite fraction af- 
fecting the results, we fix the satellite fraction to the value 
we find by fitting the halo model to all lenses in each stellar 



18 



Edo van Uitert et al.: The relation between baryons and dark matter 




Fig. 18. The halo mass as a function of the mean stel- 
lar mass for early-type galaxies in different redshift slices. 
Although we lack the statistical precision to draw definite 
conclusions, the measurements support the view that at 
the high stellar mass end, galaxies at a higher redshift have 
lower halo masses. 



mass bin. We apply the various corrections (e.g. scattering 
of lenses between bins), and show the results in Figure [T51 
The errors on the best fit halo masses are large, and 
we therefore do not obtain tight constraints on the evolu- 
tion of the halo masses for the low stellar mass bins. For 
the highest stellar mass bin, however, it appears that the 
halo mass is smaller by roughly a factor of two for the two 
highest redshift slice s. The redshi f t depe ndent stellar-to- 
halo mass relation of iMoster et al.l (J2010[ ) predicts that at 
M* = 6 X 10^^ Mq, the halo mass increases by ^ 35% be- 
tween z = 0.5 and z = 0.0. In lLeauthaud et all (|201l[ ). the 
evolution of the stellar-to-halo mass relation from z = 1 
to z = 0.2 is studied using a combined galaxy-galaxy weak 
lensing, galaxy spatial clustering, and galaxy number densi- 
ties analysis in the COSMOS survey fScoville et al., .200,7i) . 
At stellar masses M* > 10^^, the halo mass appears to 
decrease with redshift for a given stellar mass, but the 
small v olume probed by C OSMOS prevents a clear de- 
tection. I Brown et al.l (|2008f ) study the growth of the dark 



matter content of massive early-type galaxies between a 
redshift of 0.0 and 1.0 by measuring the space density 
and spatial clustering of the galaxies. They find that be- 
tween redshift z = 1.0 and z ~ 0.0, the dark matter haloes 
grow with ~100%, while t he stellar masses of these galax- 
ies only grow with ^30%. IConrov et al.l (l2007l ) utilizes the 
motions of satellite galaxies around isolated galaxies to con- 
strain the evolution of the virial-to-stellar mass ratio, and 
they find that between z ~ 1 and z ^ this ratio re- 
mains constant for host galaxies with a stellar mass below 
1.5 X 10^ ^Af©, but increases by a factor 3.3 ± 2.2 for hosts 
with M* > 1.5 X IO^^Mq. These findings are in qualitative 
agreement with our results. 



8. CONCLUSIONS 

We measured the halo masses for early- and late- type galax- 
ies and compared these to their luminosity and stellar mass. 
For this purpose, we measured the weak lensing signal in- 
duced by the galaxies with SDSS spectroscopy that overlap 
with the RCS2, and modelled the data with a halo model. 
This enabled us to improve the constraints on the lensing 
measurements for the most massive galaxies, which typi- 
cally reside at redshifts where the SDSS is not very sensi- 
tive. 

The halo mass and the dynamical mass correlate well for 
early- type galaxies, but not for late- type galaxies. A likely 
explanation is that late- type galaxies are rotating, resulting 
in an overestimation of the velocity dispersion, and hence 
of the dynamical mass. Furthermore, in contrast to the dy- 
namical mass, the weak lensing mass can easily be related 
to numerical simulations, and provides constraints for the 
models that describe the relationship between baryons and 
dark matter. 

The halo masses of galaxies increase with luminosity 
and stellar mass. For a given luminosity, the halo mass 
of the early-types is on average about five times larger 
than the late- types. We fitted a power law relation between 
the luminosity and halo mass, and find that in the range 
10^*^ < Lr < 10^^-^ Lq, the halo mass scales with lu- 

9 04+0.09 9 + 0.7 

minosity as M/j ex L *-o.i6 and M/j ex L -o^ for the 
early- and late-type galaxies respectively. For an early-type 
galaxy with a fiducal luminosity Lq = lO^^L^,©, we obtain 
a mass M200 = {1.93^0,11) x IO^^/i-^Mq. We computed 
L2001 the additional luminosity around the lenses within 
r200) and find that the M200 7^200 ratio of the early- types 
is larger than for the late-types: for ^200 < lO^^i© we find 
^200/^200 = 42 ± 10 for early-types, whilst M2oo/i20o = 
17 ± 9 for late- types. This suggests that the difference in 
halo mass is not solely due to the fact that early-types re- 
side in denser environments, but is at least partly intrinsic. 

Below a stellar mass of lO^^M© the halo mass of early- 
and late-types are comparable. For larger stellar masses, 
the best fit halo masses of the early-types are larger than 
the late-types. We computed Af*.200i the total stellar mass 
within r2oo ; in order to calculate the baryon conversion effi- 
ciency rj. Our results for early-type galaxies suggest a varia- 
tion in efficiency with a minimum of ~10% for a stellar mass 
-^^*,2oo = lO^^M©. The results for the late-type galaxies are 
not well constrained, but do suggest a larger value. 

The satellite fraction is ^40% for the low luminosity 
(stellar mass) early-type galaxies, and decreases rapidly 
to < 10% with increasing luminosity (stellar mass). The 
satellite fraction of the late-types has a value in the range 
0-15%, independent of luminosity or stellar mass. The 
satellite fraction is difficult to constrain at the high stel- 
lar mass/luminosity end, as the shape of the combined 
shear signal from the satellites mimics an NFW profile. 
Decreasing the truncation parameter leads to tighter con- 
straints, and appears to be justified for the most massive 
ea rly-type satellite s based on the N-body simulation results 
of iLimousin et al.l (|2009( ). Additional support comes from 
studying the shear signal of massive early-type galaxies that 
were selected to be satellites shown in Appendix lUl but the 
errors are currently too large to constrain the fraction of 
dark matter that is stripped. A more realistic description 
of the stripping of the haloes of massive satellite galaxies 
may result in an improvement of the constraints on the 
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satellite fraction from weak lensing studies alone. 

The halo mass appears to decrease with redshift for 
the highest stellar mass bins, a trend that is qualitatively 
in agreement with predictions from numerical simulations. 
The signal-to-noise on the measurements is currently too 
low to provide a detailed view on the growth of dark matter 
haloes, but it shows that with future surveys weak lensing 
can be used to study in great detail the evolution of the 
relation between baryons and dark matter. 
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Appendix A: Scatter of lenses between bins 

In this Appendix we describe how we calculate the bias 
that resuhs from the scatter of galaxies between lensing 
bins due to the stellar mass errors. The bias that resuhs 
from the scatter due to the luminosity errors has been 
calculated in a similar fashion. To begin, we create a large 
set of simulated lens catalogues. We construct the stellar 
mass function from the MPA/JHU catalogue, randomly 
draw stellar masses from this distribution and assign these 
to our lenses. We fit a powerlaw of the form A/200 = a*Aff * 
to our initial observations, and calculate the halo mass 
of each galaxy. Assuming that the density profile of each 
lens follows an NFW profile, we calculate the ellipticities 
of the source galaxies under the assumption that they are 
intrinsically round. Next we create 20 new lens catalogues 
by applying a log- normal scatter with a width of 0.1 to 
the stellar masses. We use the stellar mass bins from Table 
[3] to stack the lensing signal of the scattered lenses, and 
measure the tangential shear using the original source 
catalogue. We fit the lensing signal between 30 and 200 
kpc with the halo model, imposing a satellite fraction of 
0% as the lenses were randomly inserted in the images. 
The ratio of the best fit halo masses for the original lenses 
and the lenses with scattered stellar masses gives the bias. 
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Fig. A.l. The ratio of the best fit halo mass for the un- 
scattered lens catalogue to the best fit halo mass for the 
lenses to which a log-normal scatter of 0.1 in stellar mass 
has been applied, for the early- type galaxies {left) and the 
late- type galaxies (right). The halo masses are underesti- 
mated at the high stellar mass end due to low mass objects 
scattering into the high mass bins. 



As the stellar mass function and the best fit pow- 
erlaw are different for the two galaxy types, we make 
two sets of simulations to study the bias for early- and 
late-type galaxies separately. We do not account for 
evolution with redshift, although the stellar mass function 
evolves between z = 0.0 and z = 1.0, most strongly for 
Mh, < 10"Afo (e.g. IVulcani et all 120111: iPozzetti et all 
l2010l ). We are only sensitive to the change of the shape 
of the stellar mass function, which is most noticeable for 
10^°-^ < M* < lO^^M©. However, the bias in this regime 
is small, and we do not expect the change in shape to 
strongly affect our results. The relation between stellar 
mass and hal o mass may also evolve between z = 0.0 an d 
z = 0.5 (e.g. iMoster et aP I2OIOI: iLeauthaud et al.ll201lD . 
We currently lack sufficient signal-to-noise to study this 
in detail. As we will demonstrate, the bias is not very 
sensitive to changes in the powerlaw slope, and a mild 
evolution does not significantly alter the results. 

The ratio of the input halo masses to the best fit halo 
mass measured for the lenses that have been scattered is 
shown in Figure lA.ll The error bars indicate the standard 
deviation of the simulations. We find that the bias is 
highest for the early-types at the high mass end. This is 
due to the steepness of the stellar mass function, which 
leads to the net effect that low stellar mass objects scatter 
into and contaminate the high stellar mass bins. The bias 
for early-types at the low mass end is slightly smaller than 
1, as the stellar mass function turns over at ^ 5 x 10^°M* 
and becomes smaller with decreasing stellar mass. The 
stellar mass function of the late-types is monotonically 
decreasing, and consequently the bias docs not become 
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smaller than unity. At the high mass end, the stellar mass 
function of the late-types is poorly determined due to the 
lack of objects. We cannot reliably determine the bias for 
the S6 late-type bin, and therefore apply the correction 
factor of the S5 bin to this bin as well. 

The observed stellar masses have already been scat- 
tered, and the best fit powerlaw is therefore too shallow. 
To investigate how this affects the bias, we correct our 
initial halo masses for the scattering, and again fit a 
powerlaw between stellar mass and halo mass. We repeat 
our simulations with these new powerlaw slopes, and 
find that the correction factors change by at most 4%. 
The correction we apply is obtained using the corrected 
powerlaw slopes. 

The intrinsic stellar mass function is steeper than the 
observed one as on average more low stellar mass objects 
have scattered upward. Although we cannot retrieve the 
intrinsic stellar mass function, we can obtain an estimate 
of the level of contamination. For this purpose, we draw 
1 X 10® objects from the observed stellar mass function, 
apply the log-normal scatter, and compare the number of 
objects in the stellar mass bins before and after the scatter. 
The number of lenses in the three lowest stellar mass bins 
does not change much after the scatter, but it increases 
with stellar mass for the more massive bins, reaching a 
maximum of 36% more lenses in the S7 early-type bin. The 
increase in the number of objects may be even larger, as the 
observed stellar masses have already been scattered, and 
therefore the observed stellar mass function is smoother 
than the intrinsic one. As the stellar mass function at the 
high mass end is already very uncertain, we do not attempt 
to retrieve the intrinsic stellar mass function. However, 
the bias correction is sensitive to the slope at the high 
mass end, and the correction factors may actually be larger. 



Appendix B: Mean versus fitted halo mass 

The distribution of halo masses for a certain luminosity (or 
stellar mass) is given by the conditional probability func- 
tion, which is usually described by a log-normal function of 
the form 



P{mh\l) ex exp ( - 






(B.l) 



where I — log(L), m/i — \og{Mh) and ct™^ is the scatter 
in nih. In this Appendix we study how the best fit lensing 
mass is related to either the mean halo mass or to the 
centre of the halo mass distribution, nih^cent- To mimic 
the selection of real galaxies, we assign a value to mh^cent 
and CTm^ , and randomly draw 1000 galaxies from the 
conditional probability function which has been convolved 
with the halo mass function fEauation [T0|) . We calculate 
the NFW shear profiles of these galaxies, average their 
signals to simulate the usual lensing procedure, and fit an 
NFW profile to the stacked shear. Figure IB. II shows the 
ratio of Mh.cent to the best fit NFW mass in the top panel, 
and the ratio of the mean halo mass to the best fit NFW 
mass in the lower panel. The lines correspond to different 
values of fTmh, ranging from 0.10 to 0.40 from bottom to 
top. Note that the scale of the vertical axes in the two 
panels is different. 

In Figure IB. lb . we see that the best fit NFW mass is 
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Fig. B.l. The ratio of the central mass of the halo mass 
distribution, mh^cent, and the best fit NFW mass (top) and 
the ratio of the mean halo mass and the best fit NFW mass 
{bottom) as a function of best fit NFW mass. Different lines 
correspond to values of am^ 0.10 (bottom line), 0.15, 0.20, 
0.25, 0.30, 0.35 and 0.40 (top line). The lensing mass is 
converted to the mean halo mass using the corrections from 
the bottom panel. 



considerably lower than the central mass of the distribu- 
tion. This is mainly the result of the declining halo mass 
function, which leads us to preferentially pick lower mass 
haloes. The shape of an NFW profile changes with halo 
mass because the NFW concentration parameter depends 
on halo mass. The shape and amplitude of the stacked 
shear signal is therefore not equal to the profile of an NFW 
with a corresponding mean halo mass. Therefore, the best 
fit NFW mass underestimates the mean halo mass, as 
demonstrated in Figure iBTTb . 

The ratios in Figure [BJk and IB. lb are sensitiv e to th e 
value of fTm^- We use the results from iMore et al.l ()2011[ ). 
who studied the distribution of halo masses as a function 
of luminosity and stellar mass using the kinematics of 
satellite galaxies orbiting central galaxies. As only central 
galaxies are considered in their work, the actual scatter 
for a sample of galaxies consisting of both centrals and 
satellites may be larger. On the other hand, part of the 
scatter may be introduced through uncertainties in the 
determination of the halo masses, which would imply a 
lower intrinsic scatter. 

We use Figures 4 and 9 from iMore et aO (|201l[) to 
read off the values we assign to (Tm,, for the luminosity 
and stellar mass bins. We list these values, and the 
corresponding correction factor to the mean halo mass, 
in Table IB. II The luminosities and st ellar masses in our 
sample extend to higher values than iMore et al.l ()2011h 
use, but their figures suggest that amn does not change 
rapidly at the high mass/luminosity end, and we therefore 
assume that the values remain constant. For the stellar 
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Table B.l. The values of (j,n^ assigned to the luminosity 
and stellar mass bins, and the correction factors fcorr we 
apply to convert the measured lensing mass into the mean 
halo mass. 



Sample 


o-mfe (early) 


/corr (early) 


o-m^ (late) 


/corr (late) 


LI 


0.20 


1.07 


0.25 


1.10 


L2 


0.25 


1.11 


0.29 


1.14 


L3 


0.30 


1.17 


0.30 


1.15 


L4 


0.33 


1.20 


0.33 


1.20 


L5 


0.37 


1.26 


0.34 


1.21 


L6 


0.39 


1.28 


0.35 


1.23 


L7 


0.40 


1.28 


0.35 


1.22 


L8 


0.40 


1.27 


0.35 


1.23 


SI 


0.15 


1.03 


0.10 


1.01 


S2 


0.18 


1.06 


0.10 


1.02 
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0.26 


1.12 


0.10 


1.02 
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0.32 


1.19 


0.10 


1.02 
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0.36 


1.24 


0.10 
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1.28 
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masses we use the NFW masses that have been corrected 
for the scattering of objects between the bins. 

There are further sources of uncertainty to consider 
in future studies, and we list a few of them: luminosity 
bins have a certain width, the luminosity function is not 
constant inside a luminosity bin, and lens galaxies are 
located at a range of redshifts. We expect that these 
complications further broaden the conditional probability 
function, which means that the correction factors we use 
may be too low. These complications should be taken 
into account to enable a detailed comparison between 
observations and simulations. 



Appendix C: Constraints on the satellite fraction 
at high halo masses 

The satellite fraction is not well constrained at the high lu- 
minosity/stellar mass end. The reason for this is illustrated 
in Figure IC.ll In Figure IC.lb we show the lensing signal of 
the L6 luminosity bin, together with the five terms of the 
halo model, using the standard truncation radius of 0.4r2oo 
for the satellite galaxies. The satellite shear signal on scales 
< 1/i^^Mpc in the halo model is the sum of stripped satel- 
lite term and the ^Wat term. It is clear that the shape of 
the combined signal is very similar to the shape of the shear 
signal coming from the central halo. As a result the error 
on the satellite fraction is large. The satellite fraction and 
the halo mass are anti-correlated, as we can see from Figure 
IC.2I The model either prefers a large mass and small satel- 
lite fraction, or a small mass and large satellite fraction. To 
reduce any bias in the best fit halo mass, we decrease the 
allowed range for the satellite fractions to a uniform prior 
between 0% and 20% for the highest stellar mass and lumi- 
nosity bins, as almost all of the galaxies in these bins are 

expected to be cent rals. 

Recent work by [Limousin et al.l (l2009l ) shows that the 
half mass radius of a subhalo is a strongly decreasing func- 



Fig. C.l. The lensing signal of the L6 early- type bin, shown 
together with the five components of the best fit halo model. 
In the upper panel the truncation radius of the stripped 
satellites is 0.4r2oo, and the shape of the combined satellite 
1-halo terms mimicks the shape of the central NFW term. 
In the lower panel the truncation radius is 0.2r2oo, changing 
the shape of the combined satellite 1-halo terms. Note that 
the halo model in the lower panel is not a fit, but serves 
to illustrate the effect of choosing a different truncation 
radius. 



tion of projected cluster-centric distance. Furthermore, the 
radial distribution of early-type satellites is more peaked 
around the cluster cen tre than the rad ial distribution of 
late-type satellites (e.g. lAnn et al.ll2008l ). Hence we expect 
that the massive elliptical satellite galaxies, which prac- 
tically always reside close to the centre of a cluster, are 
stripped of a far larger fraction of their dark matter. 

To determine whether we can observe a change in the 
truncation radius of massive early- type satellite galaxies, we 
make a selection of galaxies that are likely to be satellites 
and study their shear profile. We consider early-types in 
the mass range 10^°-^ < M* < 10^^'^^ M©, and divide them 
in three mass bins; galaxies more massive than 10^^'^^ M© 
will almost exclusively be central galaxies and hence not 
significantly stripped. To determine whether the galaxies 
are satellites or centrals, we use the SDSS DR7 photomet- 
ric redshift catalogue Photoz, which contains the photo- 
metric redshifts of 260 million galaxies, and match them to 
our source galaxy catalogue. The lenses that have a neigh- 
bouring galaxy of the same luminosity or brighter within 
750 kpc, and lie within the la errors of the photometric 
redshift of the source, are selected for the satellite sample. 
The galaxies that do not have brighter neighbours within 
1 Mpc and within the la errors of the photometric redshift 
are selected for the central sample. Note that we do not 
aim to obtain samples that are complete, but we strive to 
make a selection that enables us to quantitatively study the 
differences in the lensing signal. 
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Fig. C.2. The x^ values of the halo model fits to the L6 
early- type bin. The green star indicates the best fit. The 
three contours show the 67.8%, 95.4% and 99.7% confidence 
intervals (Ax^ of 2.3, 6.2 and 11.8 respectively). The best 
fit halo mass is anti-correlated with the best fit satellite 
fraction. 



In Figure IC.3I we show the stacked shear signal of the 
galaxies, for the central sample and for the satellite sam- 
ple, together with their halo model components. The shear 
signals of the central sample are indeed described well by 
an NFW profile. The galaxies preferentially live in isolated 
environments, and consequently the 7^ccnt term is overes- 
timated. The halo model fits of the satellite sample are 
dominated by the satellite terms, as can be observed from 
the best fit satellite fractions indicated in the plot. The 
shear signal around 100 kpc is different from the signal of 
the central sample at the same scale, and is lower than the 
halo model fit. This suggests that additional stripping of 
dark matter occurs at small scales. The measurements are 
too noisy, however, to constrain which fraction of the dark 
matter haloes is stripped. 

To illustrate the impact the choice of truncation ra- 
dius has on the best fit satellite fractions, we also con- 
sider stripped satellite profiles with a truncation radius of 
0.2r2oo- In Figure lCTb we show the shear signal of the same 
bin, but with this smaller truncation radius. Note that the 
halo model parameters are identical in both panels for il- 
lustrative purposes, and that the model in the lower panel 
is not a fit. The shear signal of the satellites at small scales 
is now clearly different from the central halo term, and the 
satellite fraction can be better constrained. We have also 
fit halo models with a truncation radius of 0.2r2oo to the 
four most luminous early-type bins. The constraints on the 
satellite fraction for both models are shown in Figure IC.4I 
The satellite fraction is better constrained for the models 
with a truncation radius of 0.2r2oo- Setting the truncation 
radius to 0.2r2oo is a rather arbitrary choice, however, and 
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Fig. C.3. The lensing signal AE as a function of physical 
distance from the lens. The lensing signal is measured for 
central galaxies (left) and for the satellite galaxies (right), 
for the 10i°-5 < M, < 1O"-°M0 bin (top), 10"-° < M, < 
lO"-5Af0 bin {middle) and 10"'^ < M, < IO^^-'^^Mq bin 
(bottom). Indicated in each plot is the number of lenses, the 
logarithm of the best fit halo mass and the best fit satellite 
fraction. The shear signal is reduced at large lens-source 
separations for the central galaxies, indicating that they are 
isolated. At small lens-source separations the shear signal of 
the satellite sample appears to be reduced compared to the 
central sample. Note that the 2-halo terms are not shown 
for clarity. 



in future studies it is necessary to include a more realistic 
pre scription for the stripping of the satellites. 

iMandelbaum et aLl ()2006() study the environmental de- 
pendence of the shear profile as a function of luminosity. 
They distinguish galaxies residing in a high-density envi- 
ronment and in a low-density environment. The brightest 
galaxies of their low-density sample are almost exclusively 
centrals, whilst in the high-density sample they are a mix- 
ture of centrals and satellites. As the lensing signal is then 
an average of the shear profiles from satellites and centrals, 
this may explain why they do not observe a reduction of 
the signal at small scales. Note that the satellite galaxies we 
study are more massive, and hence are expected to reside 
in denser environments where the haloes are stripped of a 
larger fraction of their dark matter content. 
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Fig. C.4. The best fit satellite fraction for the four highest 
luminosity bins of early- types. The thick solid (thin dashed) 
lines indicate the results determined using a truncation ra- 
dius of 0.2r2oo (0.4r2oo)- Decreasing the truncation radius 
tightens the constraints on the satellite fraction. 
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